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ABSTRACT 

We investigate shock structure and particle acceleration in relativistic magnetized collisionless 
electron-ion shocks by means of 2.5D particle-in-cell simulations with ion-to-electron mass ratios 
(mj/m e ) ranging from 16 to 1000. We explore a range of inclination angles between the pre-shock 
magnetic field and the shock normal. In "subluminal" shocks, where relativistic particles can escape 
ahead of the shock along the magnetic field lines, ions are efficiently accelerated via a Fermi-like 
mechanism. The downstream ion spectrum consists of a relativistic Maxwellian and a high-energy 
power-law tail, which contains ~ 5% of ions and ~ 30% of ion energy. Its slope is —2.1 ± 0.1. The 
scattering is provided by short-wavelength non-resonant modes produced by Bell's instability (Bell 
2004), whose growth is seeded by the current of shock-accelerated ions that propagate ahead of the 
shock. Upstream electrons enter the shock with lower energy than ions (albeit by only a factor of 
~ 5 -C mi/m e ), so they are more strongly tied to the field. As a result, only ~ 1% of the incoming 
electrons are Fermi-accelerated at the shock before being advected downstream, where they populate 
a steep power-law tail (with slope —3.5 ±0.1). For "superluminal" shocks, where relativistic particles 
cannot outrun the shock along the field, the self-generated turbulence is not strong enough to permit 
efficient Fermi acceleration, and the ion and electron downstream spectra are consistent with thermal 
distributions. The incoming electrons are heated up to equipartition with ions, due to strong electro- 
magnetic waves emitted by the shock into the upstream. Thus, efficient electron heating (> 15% of 
the upstream ion energy) is the universal property of relativistic electron-ion shocks, but significant 
nonthermal acceleration of electrons (> 2% by number, > 10% by energy, with slope flatter than 
—2.5) is hard to achieve in magnetized flows and requires weakly magnetized shocks (magnetization 
a < 10 -3 ), where magnetic fields self-generated via the Weibel instability are stronger than the back- 
ground field. These findings place important constraints on the models of AGN jets and Gamma Ray 
Bursts that invoke particle acceleration in relativistic magnetized electron-ion shocks. 
Subject headings: acceleration of particles — cosmic rays — galaxies: jets — gamma-ray burst: general 
— shock waves 



1. INTRODUCTION 

Nonthermal emission from Pulsar Wind Nebulae 
(PWNe), jets from Active Galactic Nuclei (AGNs), 
Gamma-Ray Bursts (GRBs) and supernova remnants 
(SNRs) is usually modeled as synchrotron or inverse 
Compton radiation from a power-law population of elec- 
trons accelerated at collisionless shocks. The slope of 
the power-law tail and the acceleration efficiency, i.e., 
the fraction of particles and energy stored in the tail, 
are usually parameterized ad hoc to fit the observations, 
due to the lack of a fully self-consistent theory of particle 
acceleration in shocks. 

When encountering a collisionless shock, charged par- 
ticles may be accelerated by means of two basic mech- 
anisms. In first-order Fermi acceleration (or Diffusive 
Shock Acceleration, DSA), particles stochastically dif- 
fuse back and forth across the shock front and gain en- 
ergy by scattering from ma gnetic turbulence embedded 
in the converging flows (e.g-HBlandford fc Ostriker|1978 ' 



Bell|1978[[Drury|1983|[Blandford k Eichler|l987| ). Alter 



natively, in magnetized flows, particles may also gain en- 
ergy directly from the background motional electric field 
E = —j3 x B while they gyrate around the shock. The 
latter process is given different names depending on the 
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barrier that reflects the particles back upstream, thereby 
allowing for multiple energizations. If the reflecting bar- 
rier is magnetic, due to the shock compression of the up- 
stream field, the acceleration mechanism is called Shock 



Drift Acceleratio n (SPA; e.g., Chen fc Armstrong||1975 i 
Webb et al.|1983||Begelman fe Kirk|l990[ ). If the barrier 



is electrostatic, caused by electron-ion charge separation 
at the shock, the p rocess is named Shock Surfing Accel- 
eration (SSA; e.g., Lee et al. 1996 Hoshino & Shimada 



2002 



Shapiro & Ucer 2003 1 . The SSA mechanism can 



only operate in electron-ion shocks (and not in electron- 
positron shocks), since no electrostatic barrier appears if 
the incoming species have the same rigidity. In general, 
not only the elemental composition, but also the mag- 
netization and bulk Lorentz factor of the upstream flow 
may be important in determining the mechanism respon- 
sible for particle energization and the resulting acceler- 
ation efficiency. If the upstream medium is magnetized, 
an additional parameter is the obliquity angle that the 
upstream magnetic field makes with the shock direction 
of propagation. 

In the absence of turbulence, cold particles are con- 
strained to slide along the ordered magnetic field lines, 
which are advected downstream from the shock. For high 
magnetic inclinations, in order to return upstream the 
particles shoul d be moving along the fiel d faster than the 
speed of light (Begelman & Kirk 1990). If efficient ac- 
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celeration requires repeated crossings of the shock, such 
"superluminal" geom etries should be very inefficient for 



particle accele ration (Kirk & Heavens 1989 Ballard & 



Heavens|1991 ). This conclusion can be questioned in the 



following cases: (i) a fraction of the incoming particles 
is pre-heated such that their gyro-radius is larger than 
the thickness of the shock, and they can experience mul- 
tiple shock crossings before being advected downstream; 
(ii) there are appreciable fluctuations of the magnetic 
field at the shock, that may create local "subluminal" 
configurations where particle acceleration can occur. In 
both cases, the final outcome is difficult to predict. In 
support of the pre-heating scenario, a candidate mecha- 
nism has been proposed for electron heating in magne- 
tized oblique shocks. The s ynchrotron maser ins tability 
at relativistic shock fronts (Langdon et al. 1988) gener- 
ates a coherent train of intense electromagnetic precur- 
sor" waves propagating upstream, which push on elec- 
trons and make them lag behind ions. The difference 
in bulk velocities between electrons and ions generate s 
longitudinal electrostatic oscillations (Lyubarsky 2006), 
that e ventually dissi pate by heating the upstream elec- 
trons ( |Hoshino|2008[ ). However, it is not clear if this will 
result in electron acceleration, besides heating. More- 
over, if the incoming electrons are too hot, the syn- 
chrotron maser instability may be suppressed, thereby 
self-limiting the whole process. 

Alternatively, as mentioned above in (ii), particle ac- 
celeration in superluminal shocks may be permitted if 
a sufficiently strong turbulence reorients the field at the 
shock, opening subluminal channels where particles can 
return upstream. If the pre-shock medium is not turbu- 
lent by itself, such magnetic field fluctuations need to be 
generated by shock-accelerated particles that escape up- 
stream across field lines. However, the existence of these 
particles is not obvious a priori, since in turn it requires 
sufficient magnetic turbulence to mediate their accelera- 
tion. In any case, the highly nonlinear problem of wave 
generation and particle heating and acceleration needs 
to be addressed with a self-consistent approach. 

Fully kinetic particle-in-cell (PIC) simulations provide 
a powerful tool for exploration of the structure of col- 
lisionlcss shocks from first principles, thus determining 
self-consistently the interplay between shock-generated 
waves and accelerated particles. As opposed to semi- 
analytic kinetic theory methods (e.g., |K irk fc Heavens 

T> 



19891 [ Ballard fe Heavens 



Ostrowski 2004 



,19911 |Achterberg et al . , 
Keshet & Waxman 2005 or Monte Carlo test-particle 



simulations ( e.g. ," Bednarz fc Ostrowski|1 998; Nic miec "fc l 
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PIC simulations can tackle the problem of particle ac- 
celeration in shocks without the need for simplifying as- 
sumptions about the nature of the magnetic turbulence 
or the details of wave-particle interactions. 
Multi-dimensional PIC simulations of relativistic un- 



magne tized shocks have been presented by Spitkovsk 



( 2005 1 for electron-positron flows (see a l so Chang et a 
20UgPKeshet et al.| [2009] |Haugboelle| |2010J, and 



Spitkovsky| ( | 2008a |) for elec tron-ion flows. Spitkovsky 



by 



(20086) and Martins et al.| (|2009j) have shown, respec- 
tively for electron-positron and electron-ion plasmas, 
that unmagnetized shocks naturally produce accelerated 
particles as part of the shock evolution. A comprehensive 
study of particle acceleration in relativistic magnetized 



electron-positron sh ocks has been performed by |Sironi fc| 
Spitkovsky (2009a hereafter SS09) by means of multi- 
dimensional PIC simulations. In agreement with earlier 
one-dimensional (ID) results ( |Langdon et al.fl988 Gal 



lant et al.|1992[ ) , they found that particle acceleration m 



superluminal pair shocks is extremely inefficient, mean- 
ing that self-generated turbulence is not strong enough 
to allow for significant particle diffusion across field lines. 

In this work, we extend the analysis of SS09 by in- 
vestigating via 2.5D PIC simulations the properties of 
relativistic magnetized electron-ion shocks. In particu- 
lar, we explore how the level of electron heating (namely, 
the fraction of pre-shock bulk kinetic energy transferred 
to post-shock electrons) and the efficiency of particle ac- 
celeration depend on the upstream bulk Lorentz factor, 
magnetization, and magnetic obliquity. We find that 
electrons are heated up to a few tens of percent of the 
upstream ion energy, regardless of the upstream condi- 
tions, and they can even reach equipartition with ions in 
high-obliquity magnetized shocks. With regards to parti- 
cle acceleration, efficient nonthermal energization of ions 
(and electrons, to a lesser degree) via a Fermi-like process 
is observed for subluminal configurations. The pressure 
of shock-accelerated ions propagating ahead of the shock 
can substantially perturb the incoming flow and alter 
the shock structure. In contrast, the downstream spec- 
trum of both ions and electrons in superluminal shocks 
does not show any signature of acceleration to nonther- 
mal energies, in agreement with the results of SS09 for 
electron-positron shocks. 

This work is organized as follows. In Sj2]we discuss the 
setup of our simulations and the magnetic field geom- 
etry. In fj3 the shock structure and internal physics is 
investigated for one representative subluminal and one 
superluminal obliquity; for the same angles, in Qwe ex- 
plore the mechanisms responsible for particle heating and 
acceleration. The main results of our work, concerning 
the efficiency of electron heating and particle accelera- 
tion as a function of bulk Lorentz factor, magnetization 
and magnetic obliquity, are presented in ^5] We summa- 
rize our findings in ^6] and comment on the application 
of our results to astrophysical scenarios. 

2. SIMULATION SETUP 

W e use the 3D elec tromagnetic PIC code TPJSTAN- 
MP ( |Spitkovsky||2005 1 , which is a parallel ver sion of the 
publicly available code TRISTAN ( |Buneman|l993[ ) that 
was optimized for studying collisionless shocks. Our sim- 
ulation setup parallels SS09 very closely, which we repeat 
here for completeness. 

The shock is set up by reflecting a cold electron-ion 
"upstream" flow from a conducting wall located at x — 
(Fig. [I]). The interaction between the incoming beam 
(that propagates along —x) and the reflected beam trig- 
gers the formation of a shock, which moves away from 
the wall along +x. This setup is equivalent to the head- 
on collision of two identical plasma shells, which would 
form a forward and reverse shock and a contact discon- 
tinuity. Here, we follow only one of these shocks, and re- 
place the contact discontinuity with the conducting wall. 
The simulation is performed in the "wall" frame, where 
the "downstream" plasma behind the shock has zero x- 
velocity. Since the downstream plasma is at rest, no extra 
boosts are needed to study the spectra. 
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We perform simulations in both 2D and 3D compu- 
tational domains, and we find that most of the shock 
physics is well captured by 2D simulations. Therefore, 
to follow the shock evolution for longer times with fixed 
computational resources, we mainly utilize 2D runs. All 
three components of particle velocities and electromag- 
netic fields are tracked, however. So, our simulations are 
effectively "2.5D", i.e., 2D in physical space but 3D in 
momentum space. 

We use a rectangular simulation box in the xy plane, 
with periodic boundary conditions in the y direction 
(Fig.nl). Each computational cell is initialized with two 
electrons and two ions, but we also performed limited 
experiments with a larger number of particles per cell 
(up to 8 per species), obtaining essentially the same re- 
sults. The relativistic electron skin depth for the incom- 
ing plasma (c/cj pe ) is resolved with 10 computational 
cells and the simulation timestep is At = 0.0450;"^. 

Here, w po = {Aire 2 n e /~fom e ) 1 f 2 is the relativistic elec- 
tron plasma frequency for the upstream flow, with elec- 
tron number density n e (measured in the wall frame) and 
bulk Lorentz factor 70 . We employ a reduced mass ra- 
tio mi/m e = 16, which allows to follow the shock evolu- 
tion for longer times (in units of uj^ 1 = (mj/m e ) 1/ ' 2 Wp e 1 ), 
while still clearly separating the ion and electron dynam- 
ical scales. As we show in Appendix [XJ we obtain essen- 
tially the same results when using higner mass ratios (we 
tried m.;/m e = 100, and up to rrii/m e = 1000 in some 
cases), which suggests that a mass ratio rrii/m e = 16 is 
already "large" enough to capture the correct accelera- 
tion physics in our shocks. In units of the relativistic 
ion skin depth c/uj p i — (mi/m e ) 1 ' 2 c/w pe , the computa- 
tional domain is ~ 25c/w p j wide (along y), which corre- 
sponds to 1024 cells for rrii/m e — 16 and 3072 cells for 
mi/m e = 100. We also tried larger domains, up to 4096 
cells wide (or ~ 100c/w p i for mi/m e — 16), obtaining 
similar results. Our computational box expands along x 
at the speed of light (see below), reaching the length of 
~ 10,000 c/wpi (or - 400,000 cells for mi/m e = 16) at 
the final time of our longest simulation. 

The incoming electron-ion stream is injected along — x 
with bulk Lorentz factor 70 = 15 and a small thermal 
spread A7 = 10~ 4 . As we comment in f|5] our results 
can be rescaled to account for a different 70 (we ex- 
plore a wide range in 70, from 3 to 50). The upstream 
flow is seeded with a background magnetic field Bo such 
that the ratio of magnetic to kinetic energy density is 
a = BQ/AirjomiUiC 2 — 0.1, where rii (—n e ) is the number 
density of incoming ionsj^j In fj5]we show how our results 
change for lower magnetizations, down to a = 10~ 5 , at 
which point the shock properties approach the unmagne- 



tized case discussed by Spitkovsky (2008a). The magne 



tized regime we explore here may be relevant for internal 
shocks in GRBs and AGN jets, as we discuss in Ej6j 
For each value of the magnetization, we explore a range 

1 The magnetization parameter as defined above pertains to ions. 
For electrons, <r e o = (mi/m e )cr at injection, where the electron 
kinetic energy is 7om e e 2 . So, it seems that the electron magne- 
tization would depend on m;/m e , for fixed cr. However, on their 
way to the shock, electrons increase their average energy up to a 
fraction a of the initial ion energy ■yorriiC 2 , where a is independent 
olm,i/m e (see Appendix [A|| . So, the "effective" electron magneti- 
zation <r e c ff = (1/ct) a does not depend on the mass ratio. 




Fig. 1. — Simulation geometry. The magnetic field (red arrow) 
can be either out of the simulation plane (as sketched here) or in 
the plane. 



of magnetic inclinations by varying the angle 9 between 
the shock direction of propagation +x and the upstream 
magnetic field Bo (Fig. [T]). The magnetic obliquity angle 
9 is measured in the walT frame. We vary 9 from 9 = 0°, 
which corresponds to a "parallel" shock, with magnetic 
field aligned with the shock normal, up to 9 = 90°, i.e., 
a "perpendicular" shock, with magnetic field along the 
shock front. For 9 ^ 0°, in the upstream medium we also 
initialize a motional electric field Eo = — f3 x Bo, where 
f3 = —fio x is the three-velocity of the injected plasma. 
As we show in Appendix[B] our results do not depend on 
the orientation of the magnetic field with respect to the 
simulation plane, as parameterized by the azimuthal an- 
gle <p (cp = 0° for field in the xy plane of the simulations, 
(f = 90° for field in the xz plane). 

As a byproduct of the shock evolution, particles and 
electromagnetic waves may propagate upstream from the 
shock at the speed of light. To ensure that the long-term 
evolution of the shock is captured correctly, they should 
not be removed from the computational domain, oth- 
erwise we would artificially suppress any feedback they 
may have on the shock. For this reason, we employ a 
"moving injector" (receding from the wall at the speed 
of light) and an expanding simulation box (see SS09 for 
details) . This permits us to follow the shock evolution as 
far as the computational resources allow, preserving all 
the particles and waves generated by the shock. 

However, since the shock velocity is smaller than the 
speed of light, the distance between the injector and the 
shock increases linearly with time. In all PIC codes, 
a numerical heating instability arises when cold rela- 
tivistic plasma propagates for large dis tances over the 
numerical grid ( jDieckmann et ah 2006 ) . For electron- 
positron shocks, in order to suppress this instability we 
let our moving injector periodically jump backward (i.e., 
towards the wall), so that the distance between the shock 
and the injector stayed roughly constant (see the "jump- 
ing injector" technique in SS09). However, for electron- 
ion shocks this is not a viable solution, since the incoming 
flow is heavily influenced by particles and waves gener- 
ated by the shock, as we show in f|3] By resetting the 
position of the injector, we would discard such particles 
and waves, thus dramatically interfering with the long- 
term evolution of the shock. Instead, we notice that the 
onset of the numerical instability is significantly delayed 
if we reduce the transverse size of our simulations. For 
long (> 2000 Wpj 1 ) simulations of strongly magnetized 
high-obliquity shocks, which are more affected by the 
numerical instability, we then employ a computational 
box with only 256 transverse cells (down to 64, in some 
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cases). At Lo p it < 2000 (i.e., before the growth of the 
instability), the results from such small simulations are 
identical to our "fiducial" runs with 1024 transverse cells, 
showing that the relevant 2D properties of the shock are 
still preserved. For mj/m e = 16, a transition to the ID 
regime is observed only for computational boxes with less 
than 32 transverse cells. 

3. SHOCK STRUCTURE 

We now describe the structure and internal physics of 
relativistic magnetized electron-ion shocks. In this sec- 
tion, as well as in <HJ we fix the upstream bulk Lorentz 
factor (70 = 15) andrnagnetization (a ~ 0.1), and we in- 
vestigate two representative magnetic inclination angles: 
9 = 15° and 9 = 75°. The dependence of our findings 
on the bulk Lorentz factor, magnetization and magnetic 
obliquity of the upstream flow is discussed in Jj5] 

According to the criterion anticipated in |lj a shock 
is "superluminal" if particles cannot escape ahead of 
the shock by sliding along the magnetic field. In the 
upstream frame, this corresponds to magnetic inclina- 
tions 9' such that cos 9' > cos9' crit — (3' sh , where f3' sh 
is the shock speed in the upstream frame. In the sim- 
ulation frame, the critical obliquity angle that sepa- 
rates subluminal and superluminal configurations will be 
9 clit = arccot[7 s h (A> + Ash)]> where /3 sh and 7 sh arc the 
shock velocity and Lorentz factor in the simulation frame 
(see SS09 for a detailed discussion). For 70 = 15 and 
a = 0.1, the critical obliquity angle is 6* cr i t ~ 34°, and it 
stays confined within a relatively narrow range (between 
~ 26° and ~ 42°) for relativistic (70 > 2) flows with 
moderate magnetization (er < 1.0), as shown in Fig. 2 of 
SS09. It follows that 9 = 15° is a subluminal shock, and 
9 = 75° is superluminal. 

Figs. 2][6 present the internal structure of the shock 
as a function of the longitudinal coordinate x. For the 
subluminal angle 9 = 15°, Fig. [2] covers the whole longi- 
tudinal extent of the simulation domain, whereas Fig. [3] 
focuses on a smaller region around the shock, as delim- 
ited by the vertical dashed red lines in the first panel of 
Fig.ra Bothfi gures refer to u pl t = 2250, when the shock 
is already fully developed. The temporal evolution of the 
shock, from the early stages up to the self-similar state 
approached at ui pi t > 2000, is shown in Figj4j Figs. [5] 
and [6] correspond respectively to Figs. [2]and[3[ior the su- 
perluminal angle 9 — 75°, at time ui p it = 1350. Here, we 
choose an earlier time (with respect to co p it = 2250 for 
9 = 15°), to avoid the numerical instability discussed in 

§ which grows at w p it ~ 2000 in high-obliquity shocks, 
t, time u pi t — 1350 already captures the main proper- 
ties of the long-term evolution of the sho ck. 

This section is organized as follows. In §3.1| we present 
a general description of the structure of relativistic mag- 
netized shocks in electron-ion plasmas. We focus on the 
properties they share with electron-positron shocks, and 
emphasize the differences. The n, we separately discuss 
the cases of subl uminal (in S3.2 for 9 = 15°) and super- 
luminal (in S3.3 for 9 = 75°) shocks, describing how the 
obliquity of the field affects the shock structure. 

3.1. Magnetized Relativistic Electron-Ion Shocks: 
General Overview 

For all oblique magnetic configurations (i.e, with the 
exception of strictly parallel and perpendicular shocks), a 



fluid structure with two shocks (a strong "fast" shock and 
a weak "slow" shock) is seen in our PIC simulations of 
both electron-positron (SS09) and electron-ion flows, as 
expected from analytic theory ([Majorana" fc Anile|p87| 
and MHD calculations (Komissarov 2003). The incom- 
ing fluid does not stop completely at the fast shock. The 
residual bulk velocity is larger than the local slow mag- 
netosonic speed, and a slow shock is formed behind the 
fast shock (i.e., closer to the wall, in our simulations). 
The slow shock transition in our PIC simulations is not 
as sharp as MHD would predict, due to insufficient dissi- 
pation along the magnetic field. The slow shock becomes 
slower and weaker as 9 approaches either 0° or 90°, and 
relatively stronger and faster for intermediate obliquities. 
For oblique configurations, we will neglect the slow shock 
and refer to the fast shock simply as "the shock" . 

The jump in density and electromagnetic fields at the 
fast shock, as well as the shock velocity, are in agree- 
ment with MHD calculations, and with PIC simulations 
of electron-positron shocks (see Table 1 in SS09). As ob- 
served in pair shocks, at low obliquities (9 < 45°) the 
effective adiabatic index of the downstream plasma is 
4/3, as in a 3D relativistic gas, whereas for 9 > 45° it 
tends to 3/2, as in a 2D fluid. In fact, right behind the 
shock the particle motion is mostly confined to the plane 
orthogonal to the field, which for large inclination angles 
is nearly degenerate with the direction of propagation 
of the incoming flow. This prevents efficient isotropiza- 
tion along the field and results in a 2D adiabatic indexj^] 
Isotropization proceeds with distance behind the shock, 
and the particle distribution far downstream is roughly 
isotropic. 

The mechanism that mediates randomization of parti- 
cles at the shock varies depending on the field obliquity. 
In electron-positron magnetized plasmas, SS09 found 
that low-obliquity shock s are mediated by Weibel-like fil- 
amentation instabilities (Wcibcl 1959; |Medvedev fc Loeb 

Tha 



1999 Gruzinov fc Waxman|1999 ), as h appens tor unmag- 



netized shocks 



e.g. 



Spitkovsky 



2005| )p| The free energy 



for the instability comes from the counter-streaming be- 
tween the incoming flow and the shock-accelerated par- 
ticles propagating upstream. 

Low-obliquity shocks in electron-ion magnetized plas- 
mas are still mediated by counter-streaming instabilities, 
but the nature of the instability is different. For our 
fiducial values 70 = 15 and a = 0.1, the shock evolu- 
tion at early times (w p ji < 300) is governed by the elec- 
tron Weibel instability, which has the largest growth rate 



(~ uj pc ^> w P i). Later on, the so-called Bell's current- 
driven instability de velops ( |Bell fc Lucekl2001[|Beil|2004| 
Reville et al. [2006), and the shock results from nonlm- 



ear steepening of the circularly-polarized A lfven ic-type 
waves associated with the instability (see ^3.2 for de- 
tails). The source for the instability is the electric cur- 

2 We remark that this effect is not an artificial consequence of the 
reduced dimensionality of our computational domain, but it holds 
also for 3D simulations. Yet, it is more severe for 2D simulations 
with out-of-plane fields, since in this case the plane perpendicular 
to the field almost coincides with the simulation plane, and particle 
isotropization is even less efficient (see Appendix [B| . 

3 More precisely, the instability results from the coupling of the 
electromagnetic filamcntation (Weibel) instabil ity with the elec- 
trostatic two-stream instability, as explained by Bret (2009J. The 
resulting modes are oblique with respect to the streaming direction. 
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rent of shock-accelerated ions that propagate ahead of 
the shock, and the polarization of the resulting modes is 
non-resonant with respect to the ion gyro-motion. As 
we argue in Appendix O Bell's instability governs the 
evolution of the shock only for a limited range of magne- 
tizations (1(T 2 < a < 0.3, for fixed 7o = 15). If a > 0.3, 
the polarization of the dominant mode changes, and the 
waves are now generated via gyro- frequenc y resonance 
with the high -energy ions heading upstream (Kulsrud & 



45 = 90° 



Pearce 



1969| ). For a < 
namically unimportant , 

the ion Weibel instab ility (as in unmagnetized shocks, 



10 z , the magnetic field is dy- 
and the shock is mediated by 



Spitkovsky||2008a| ). 



see 

In contrast, high-obliquity shocks are mediated by 
magnetic reflection of the incoming flow off the shock- 
compressed magnetic field ( |Alsop fc Arons 1988| SS09), 
for both electron-positron and electron-ion plasmas. 
The coherent gyration of the incoming particles at the 
shock triggers the so-called sy nchrotron maser insta- 
bility (Hoshino & Arons 1991), where the ring of gy- 
rating particles in the shock transition breaks up into 
bunches of charge. The bunches radiate a coherent 
train of transverse electromagnetic waves propagating 



into t he upstream ([ G allant et al., 
1992j |Spitkovsky||2005| ), with wave 



1992 
engtn 



Hoshino et al. 
comparable to 

the Larmor radius ot incoming electrons in the shock- 
compressed fields. For electron-positron shocks, this 
"precursor wave" does not have a significant impact 
on the upstream plasma (SS09). In electron-ion flows, 
the guiding-center velocity of the incoming electrons de- 
creases, since they experience relativistic transverse os- 
cillations in the strong field of the precursor wave. Due to 
their high mass, ions are less affected by the wave, and 
proceed at close to their initial velocity. The resulting 
difference in bulk velocity between t he two species wil l 
generate a longitudinal electric field ( Lyubarsky| |2006 ) , 
so that electrons are boosted toward the shock whereas 
ions are deboostedQ As a result, energy equipartition 
between electrons and ions may be achieved even before 
the flow arrives at the shock front, as we discuss in |3.3| 
On the other hand, if the efficiency of ion-to-electron 
energy transfer in the upstream is moderate or negligi- 
ble, the incoming ions will enter the shock with bulk 
energy much larger than electrons, so that they will pen- 
etrate deeper into the shock-compressed field (roughly, in 
proportion to the ratio of ion to electron Larmor radii). 
The resulting charge separation establishes a net electro- 
static field at the shock pointing toward the upstream. 
The associated cross-shock electric potential attracts the 
incoming electrons into the downstream and decelerates 
the incoming ions, so that electrons may still end up with 
a significant fraction of the ion energy. 

3.2. Subluminal Shocks: 0° < 9 < 9 CI i t 

In this section, we describe the internal structure of 
subluminal magnetized electron-ion shocks. We take 
9 = 15° as a representative case, but our results will 
apply to the whole range of subluminal angles, unless 

4 A clarification is required here regarding our terminology: by 
"boost" or "deboost" , we mean changes in the bulk ^-momentum of 
the fluid; "heating" is used when the comoving particle distribution 
broadens; finally, "acceleration" is something pertaining only to a 
subsample of the particles, extracted from the bulk of the fluid and 
energized up to suprathermal velocities. 
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Fig. 2. — Internal structure of a 9 = 15° subluminal shock at 
u) p it = 2250. As a function of the longitudinal coordinate x, the 
following quantities are plotted: (a) ^/-averaged particle number 
density, in units of the upstream value; (b) y-averaged magnetic 
energy fraction eb = B 1 / SivyominiC 2 ; (c) y-averaged longitudinal 
electric field E x , normalized to the upstream magnetic field Bo; (d) 
mean particle energy (red for ions, blue for electrons) for particles 
moving toward the shock, in units of the bulk energy of injected 
ions (orange line includes also the reflected ions; the equivalent 
for electrons overlaps with the blue line); longitudinal phase space 
of ions (e) and electrons (f); (g)-(i) particle energy spectra (red 
for ions, blue for electrons), at three locations across the flow, as 
marked by arrows at the bottom of panel (f). 



otherwise noted. In both electron-positron and electron- 
ion plasmas, subluminal shocks are characterized by the 
presence of a diffuse stream of shock-reflected parti- 
cles that propagate ahead of the shock. The counter- 
streaming between these "returning" particles and the 
incoming flow triggers the generation of waves in the up- 
stream medium, which in turn affect the structure of the 
shock and the process of particle acceleration. 

The longitudinal ion phase space in Fig. W[e) shows 
the injected ions as a cold dense beam with 7^p X i ~ —15. 
In the shock transition layer (located at ~ 730c/w p i for 
ujpit = 2250), the incoming ions are isotropized and ther- 
malizcd. A diffuse population of high-energy ions moves 
ahead of the shock with 7i/3 X i > 0, f ollowi ng the up- 
stream magnetic field; as we show in §4.1. 2| these "re- 
turning" ions have been accelerated at the shock. In 
the upstream ion spectrum (red lines in Fig. [2jh) and 
(i)), they appear as a high-energy broad bump, whereas 
the low-energy narrow peak is populated by the incom- 
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Fig. 3. — Internal structure of a 9 = 15° subluminal shock at 
uipit = 2250, zooming in on a region around the shock, as delimited 
by the vertical dashed red lines in Fig.pfa). The shock is located 
at x ~ 730c/oj p ;. As a function of theTongitudinal coordinate x, 
the following quantities are plotted: particle number density, in 
units of the upstream value (2D plot in the xy simulation plane, 
panel (a); ID y-averaged profile, panel (b)); magnetic energy frac- 
tion <e b = B 2 /87r7om i n i c 2 (2D in (c), ID in (d)); 2D plots of B y 
(e) and B z (f), normalized to the upstream magnetic field Bo; (g) 
corresponding ID profiles, for B y (solid) and B z (dashed); (h) lD 
profile of the longitudinal electric field B x , in units of Bo; (i) mean 
particle ^-momentum (red for ions, blue for electrons) for particles 
in the bulk (i.e., neglecting accelerated particles), in units of the 
momentum of injected ions (orange line includes also the acceler- 
ated ions; the equivalent for electrons overlaps with the blue line); 
longitudinal phase space of ions (j) and electrons (k). 



ing particles. If transmitted downstream, the shock- 
accelerated ions will populate a power-law nonthermal 
tail, which is seen in the energy spectrum of Fig. [2|g) 
(red line) beyond the thermal distribution. 

In contrast, very few electrons propagate back up- 
stream from the shock (see the longitudinal phase space 
of electrons in Fig. f2^f )) . This anticipates that in sublu- 
minal magnetized shocks electrons are accelerated with 
lower efficiency than ionsj^j Despite the lack of returning 

5 The very few high-energy electrons seen with 7 e /3 xe > ahead 



electrons, the electron spectrum downstream from the 
shock (blue line in panel (g)) deviates from a purely ther- 
mal distribution, which suggests that some electrons may 
be accelerated to suprathermal energies as they cross the 
shock from upstream to downstream. We refer to §4.1| 
for a detailed discussion of electron acceleration. 

The fact that returning ions greatly outnumber elec- 
trons has multiple effects on the structure of the shock, 
which were absent in electron-positron flows. First, the 
positively-charged cloud of returning ions perturbs the 
incoming flow, causing significant transfer of energy from 
ions to electrons. At the leading edge of the population 
of returning ions (x ~ 1800c/o;pi), the incoming parti- 
cles cross a kind of double layer, where ions are boosted 
toward the shock and electrons are deboosted (see the av- 
erage particle energy in Fig.j2^d) , at x ~ 1700c/u;pi; red 
for ions, blue for electrons). Beyond that point, electrons 
are significantly heated (see Fig. [2^f ) at 800c/w p ; < x < 
1500c/w p i), possibly via small-scale plasma oscillations 
(Spitkovsky et al., in prep). The average electron energy 
increases toward the shock at the expense of the ion bulk 
energy, reaching up to ~ 15% of the kinetic energy of in- 
jected ions (blue line in Fig. [2](d)) . However, electrons 
still enter the shock with lower average energy than ions. 
The resulting charge separation in the shock-compressed 
magnetic field generates a net cross-shock electrostatic 
field (E x > 0, seen at x ~ 730 c/uj pi in Fig. gc)). The 
associated cross-shock potential further increases the av- 
erage electron energy; in the downstream, electrons have 
~ 30% of the initial ion energy. 

Furthermore, the imbalance between returning ions 
and electrons establishes a net (ion-driven) current in 
the upstream medium. The compensating current by the 
incoming electrons triggers the growth of the so-called 
Bell's instability, which is expected in this regime o f 
shock parameters (Bret 2009 Lemoine fc Pelletier|2 010). 
As shown in panel (e) (for B y ) and (t) (for B z ) of Fig.J3j 
the instability generates circularly-polarized Alfvemc- 
type waves in the upstream medium, with wavevector 
along the background fieldj^] The polarization of the 
waves is such that they are non-resonant with respect 
to the returning ions, and their wavelength (~ 25c/w p j) 
is much smaller than the Larmor radius of the highest- 
energy particles. As shown in Fig. [3^g) (solid for B y , 
dotted for B z ), the amplitude of the waves increases as 
they approach the shock, reaching up to ~ 4 times the 
strength of the initial field (measured in the simulation 
frame). Right before the shock (x ~ 730c/cj pi ), more 
than 50% of the bulk kinetic energy of the upstream 
flow has been converted into wave magnetic energy (see 
the magnetic energy fraction e B = B 2 / 8irj miriiC 2 in 
Fig. [3](d)) . After the shock, the flow isotropizes, the 
electric current required to drive the instability vanishes, 
and the waves dissipate their energy into particle energy. 
The importance of Bell's waves for ion acceleration will 

of the shock (see Fig. |2[f), and the high-energy bump in the blue 
line of Fig.|2[i)) have Deen accelerated in the first stages of evolu- 



tion, when fne shock was mediated by the electron Weibel instabil- 
ity. When ion-driven instabilities take over, electron acceleration 
shuts off. 

6 This is most clear in simulations with in-plane Bo (i.e., 
if = 0°). Here, tp = 90°, and one can only appreciate that the 
projection of the wavevector onto the simulation plane lies along x 
(see Figs.|3|e) and (f)). 
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be discussed in §4.1. 2| 

In the vicinity of the shock (730c/cj p i < x < 
790c/w p i), the pressure of returning ions is so large that 
they heavily impact the structure of the transition re- 
gion and the dynamics of the incoming flow. Due to the 
push of the returning ions, the incoming plasma slows 
down significantly, by more than 50%, well before the 
shock front (see the jump in the average x-momentum 
of incoming ions at x ~ 790c/cj p i, solid red line in 
Fig. |3l)). The deceleration and consequent compres- 
sion of the incoming fluid produces spikes in the num- 
ber density and the magnetic energy, as observed at 
x ~ 790c/o;pi in Fig. |3[a)-(d). Such spikes are tran- 
sient quasi-periodic structures, which get advected down- 
stream after ~ 100 uj p ^ (or ~ 30 w^ 1 , with the relativistic 
ion Larmor frequency ui c i = •y/o'Wpi), and then reform at 
a slightly larger distance from the shockF] This is the 
characteristic signature of the so-called "shock reforma- 
tion" process. Evidence of shock modification by the 
accelerated particles may a lso be seen in the ion energy 
spectrum, as we discuss in j 4.1.1 



(7=0.1 = 15° ^=90° 



3.2.1. Time Evolution 

The picture outlined above refers to uj p {t = 2250, when 
the number of shock-accelerated ions ahead of the shock 
is already large enough to affect the structure of the tran- 
sition region. Fig. [4] follows the shock evolution from 
u>p[t — 562 to cu pi t — 2250, showing how the y-averaged 
profiles of number density (panel (a)), magnetic energy 
fraction e B = B 2 /^TT^^rriiniC 2 (panel (b)), and transverse 
magnetic field B z (panel (c)) change with time. 

At early times (uj pl t = 562, blue line), when the pres- 
sure of returning ions is negligible with respect to the ram 
pressure of the incoming flow, the density jump occurs 
on the scale of a few ion Larmor radii, like an ideal MHD 
shock. The magnetic field at the shock (e B ~ 0.15, blue 
line in panel (b)) mostly results from the compression 
of the upstream field, with a minor contribution from 
counter-streaming instabilities. 

Starting from w p it = 1125 (yellow line), the shock 
structure changes. For 9 = 15°, the leading edge of 
the population of shock-accelerated ions recedes from the 
shock with a;- velocity ~ 0.86 c (corresponding to particles 
moving along the background field at the speed of light). 
Since the shock speed is ~ 0.32 c, the distance between 
the shock and the head of returning ions increases with 
time. Bell's instability, which at earlier times was sup- 
pressed by advection into the shock, can now amplify the 
field to highly nonlinear values (B z up to a factor of ~ 4 
larger than Bq, see green and red lines in panel (c); sim- 
ilar values for B y ). The magnetic energy at the shock 
saturates at equipartition with the upstream kinetic en- 
ergy (e B ~ 0.5, green and red lines in panel (b)). 

The increase in the strength of Bell's waves is driven 
by the growth of the population of returning ions. At 
late times (green for uj pi t = 1687, red for oj pi t — 2250), 
the shock transition region becomes much wider (panel 
(a)), due to the increasing push of the returning ions 
on the incoming flow. The time-averaged density pro- 
file (see the thin black line in panel (a) for Lo pl t = 2250) 

7 They resemble the so-called Short Large-Amplitude Magnetic 
Structures, or SLAM S, observed at the Earth's bow shock (e.g., 
|Schwartz et al.|1992[ |. 
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Fig. 4. — Time evolution of the internal structure of a 8 = 15° 
subluminal shock: ui p it = 562 (blue), ui pl t = 1125 (yellow), Lo p {t = 
1687 (green), and u) p it = 2250 (red), (a) ^/-averaged profiles of 
the particle number density; for each time, the superimposed thin 
black line is obtained by averaging over a narrow temporal window 
(lOOcjpj wide), centered on that time, (b) j/-averaged profiles 
of the magnetic energy fraction es, with the horizontal black line 
showing the value at injection (= 0.5 a = 0.05). (c) {/-averaged 
profiles of the transverse magnetic field B z (in units of Bo), with 
the horizontal black line showing the value at injection (= sin 15° ~ 
0.26). For each time, we only plot a limited region (~ 300c/w p i 
wide) around the instantaneous location of the shock. 



displays the typical structure of "cosm ic-ray modified" 
shocks (e.g., [Berezhko fc Ellison 19991, with a smooth 
density increase (the "cosmic-ray precursor") ahead of 
the main shock (here properly called "subshock"), now 
located at x ~ 730c/w p i. The magnetic energy stays 
large (e B ~ 0.5) across the whole transition region, with 
two characteristic peaks, one at the head of the cosmic- 
ray precursor, where the incoming plasma is first decel- 
erated and compressed, and the other at the subshock, 
where the flow comes to rest (see green and red lines in 
panel (b)). As clarified by panel (c), most of the mag- 
netic energy in the region between the peaks results from 
the compression of Bell's waves. 
At later times (we follow the shock evolution up to 



Upit = 5062), the shock structure does not qualitatively 
change, but the thickness of the transition region (where 
e B ~ 0.5) still expands at ~ 0.1c. A detailed study of 
the relation between shock width and mean free path of 
diffusively-accelerated ions will be presented elsewhere. 



3.3. Superluminal Shocks: 



7 crit 



< e < 90° 



We now turn to the structure of a shock with obliquity 
9 = 75°, as a representative example of superluminal 
shocks. In such shocks, particles following the magnetic 
field cannot propagate back from the shock into the up- 
stream. In fact, in Fig. [5] no particles are seen ahead of 
the shock with positive x-momentum (neither ions nor 
electrons, phase spaces in Fig.JSFe) and (f) respectively). 
The shock transition region (located at x ~ 600c/w p i 
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75° ? = 0° c5,t=1350 



Fig. 5. — Internal structure of a 



75° superluminal shock at 



ujpit = 1350. See the caption of Fig. [2]for details. Here, the slow 
shock can be seen at x < 50c/cj p i as an increase in density (panel 
(a)) and decrease in magnetic energy (panel (b)). 



for Upit = 1350), which is not perturbed by the pressure 
of energetic particles, is much thinner than for 8 = 15°, 
on the scale of a few Larmor radii of the incoming ions 
(compare the density profile in panel (b) between Fig. [3] 
and Fig. |6|. Despite the absence of returning particles, 
the incoming flow can still be perturbed by the presence 
of the shock, e.g., via electromagnetic waves emitted by 
the shock into the upst ream. 

As discussed in §3.1[ coherent Larmor gyration of the 
incoming electrons in the shock-com pressed field trig- 



gers t he synchrotron maser instability ( Hoshino & Arons 
1991 1, which creates a train of transverse electromagnetic 



waves propagating into the upstream. These "precursor" 
waves can be barely seen as short-scale ripples stretched 
alongw in the upstream region of Fig. [6tc) (for e B ) and 
Fig. Rjte) (for B y ), especially in the vicinity of the shock. 

The incoming electrons oscillate in the transverse field 
of the precursor waves, and their guiding-center velocity 
decreases. Since ions have larger inertia, their speed is 
not appreciably altered by the waves. Behind the head 
of the precursor (x ~ 1200c/w p i in Fig.j5l), the difference 
between the bulk velocities of ions an delectrons gener - 
ates a longitudinal "wakefield" i? x > ( jLyubarsky 2006 ) , 
so that electrons are initially boosted toward the shock, 
whereas ions are deboosted (see the average particle en- 
ergy in Fig. ^d) at x > 1150c/w p i; red for ions, blue for 
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Fig. 6. — Internal structure of a 6 = 75° superluminal shock at 
w pi* = 1350, zooming in on a region around the shock, as delimited 
by the vertical dashed red lines in Fig. [5|a). See the caption of 
Fig. [3] for details. 



electrons). This initiates electrostatic wakefield oscilla- 
tions in the incoming plasma (see E K wiggles in Fig. [5jc) 
at 600c/oj p i < x < 900c/w p i), that mediate a quasi- 
periodic exchange of energy between ions and electrons 
in the upstream (see analogous wiggles in Fig. [BJci) ) . 

If the strength of the precursor were uniform through- 
out the upstream region, after each wakefield oscillation 
the electron Lor entz factor would come back to the value 
70 at injection (iLyubarsky 2006). Instead, we observe a 
secular increase in the average energy of electrons (and 
decrease for ions) as the incoming flow approaches the 



shock (Fig. ffld) at 600c/w pi 



< 



< 



1100c/cj pi ). This 



corresponds To a gradient in the amplitude of the precur- 
sor waves, which are stronger in the vicinity of the shock. 
Their radiative push on the incoming electrons will then 
be larger right in front of the shock, and weaker farther 
upstream. This imbalance establishes a net wakefield 
E x > throughout the upstream region (in addition to 
the oscillations discussed above), which results in a sys- 
tematic increase of the average electron energy toward 
the shock, as seen in Fig. [5jd) . Alternatively, one could 
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define a "ponderomotive" electric potential (proportional 
to the radiative push of the precursor), whose increase 
toward the shock would ex plain the profi le of average 
electron energy in Fig.gd) ( |Hoshino[[2008| . 

As a result, electrons and ions enter the shock with 
roughly the same energy. Given the resulting absence 
of electron-ion charge separation at the shock, no sig- 
nificant cross-shock electric field is generated, so that 
the electron and ion average energies do not apprecia- 
bly change across the shock (Fig. [5f d) at x ~ 600c/u; p i). 
So, in high-obliquity shocks the upstream ponderomotive 
potential turns out to be more important than the cross- 
shock potential, yet it is not includ ed in most models o f 
relativistic magnetized shocks (e.g., Gedalin et al.|2008 ). 

The increase in electron bulk energy toward the shock 
is accompanied by substantial electron heating, i.e., 
broadening of the electron distribution. In the upstream, 
the electron spectrum shows two components (see blue 
line in Fig.[5jh) and (i)). The low-energy "thermal" part, 
that contains most of the particles, broadens and shifts to 
higher energies as the flow propagates toward the shock 
(the shift corresponds to the secular increase in electron 
energy of Fig. p[d)). The high-energy component, which 
gets more populated in the vicinity of the shock, extends 
a factor of ~ 10 higher in energy than the thermal peak. 
It is seen in the electron phase space of Fig.[5jf) as a dif- 
fuse beam moving toward the shock, with much higher 
x-momentum than the bulk flow. In the downstream 
spectrum, this extra component survives as a tail at high 
energies (blue line in Fig. | 5^ g)). 

As we discuss in §4. 2 .lj this high-energy component 
should be interpreted as a separate hotter (but still ther- 
mal) population, rather than as a nonthermal tail of ac- 
celerated particles. It is populated by electrons that, 
on their way to the shock, happen to receive a substan- 
tial kick toward the upstream by the wakefield oscilla- 
tions mentioned above. At first, the energy of these 
electrons will decrease, because the kick i s opposite to 
their i-momentum 



As explained in 4.2.2 



the electron 

energy will then oscillate with the gyro-period, periodi- 
cally reaching a maximum that, for strong kicks, may be 
much larger than the characteristic energy of electrons 
in the bulk (see Fig. [sVf ) ) . In support of this interpre- 
tation, we see in the phase space of Fig. [5jf ) that injec- 
tion of electrons into the high-energy component starts at 
x ~ 1050c/u;pi, right behind the location where the aver- 
age electron energy is driven to a minimum (Fig. p[d) at 
x ~ llOOc/oipj) by the wakefield oscillations in Fig.lMc). 
Closer to the shock, where the wakefield oscillations are 
stronger (Fig. p3Fc) at 600c/a;pi < x ~ 900c/w p i), more 
electrons are kicked back upstream and injected into the 
high-energy spectral component. This explains why the 
electron high-energy tail gets more populated toward the 
shock (compare Fig.[5jh) and (i)). In the following, since 
the high-energy component in the electron spectrum is 
ultimately powered by the dissipation of wakefield oscil- 
lations, we shall refer to this process as "wakefield heat- 
ing". 

In contrast, ions are not significantly affected by the 
wakefield oscillations, due to their larger mass. The up- 
stream ion spectrum is then consistent with a single com- 
ponent (red line in Fig. pi h) and (i)), which results in a 
purely Maxwellian distribution in the downstream region 
(red line in Fig. |5jg)). 
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Fig. 7. — Time evolution of the downstream energy spectrum in 
a 8 = 15° subluminal shock, for ions (upper panel) and electrons 
(lower panel): aj pi i = 562 (black), u pi t = 1687 (blue), ui pi t = 2812 
(green), ui p it = 3937 (yellow), and u> p it = 5062 (red). Subpanels: 
(a) time evolution of the downstream mean particle energy (red for 
ions, blue for electrons), in units of the bulk energy of injected ions, 
with the horizontal dashed blue line showing the value expected for 
electrons in the absence of any ion-to-electron energy transfer (= 
m e /mi ~ 0.06); (b) time evolution of the downstream maximum 
particle energy (red for ions; blue for electrons, multiplied by 5 for 
clarity); (c) fit to the electron (blue) and ion (red) spectrum at 
uipit = 5062, with a low-energy Maxwellian (dashed) plus a high- 
energy power law (dot-dashed). In panel (c), axes are the same as 
in the main plot. 



We remark that the main features discussed here for 
9 = 75°, and in particular the emergence of a separate 
component of hot electrons in the upstream, are common 
to all superluminal configurations, up to 9 = 90°. 

4. PARTICLE ACCELERATION 

In this section, we discuss particle acceleration in sub- 
luminal and superluminal shocks. We fix the magnetiza- 
tion (a = 0.1) and bulk Lorentz factor (70 = 15) of the 
upstream flow, and we inv estigate one representative sub- 
luminal (6 = 15°, in f 4.1 ) and one superluminal (9 = 75°, 
in [4.2) magnetic obliquity. For each case, we describe 



the time evolution of the downstream energy spectrum 
and the mechanism responsible for particle energization. 

4.1. Subluminal Shocks: 0° < 9 < 9 cr - lt 
4.1.1. Time Evolution 

Fig.[7]shows the time evolution of the ion (upper panel) 
and electron (lower panel) energy spectrum for 9 = 15°, 
in a downstream slab at fixed distance from the shock. 
As shown in panel (c), the spectrum of both ions (red) 
and electrons (blue) clearly deviates from a Maxwellian 
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(dashed), with a substantial population of nonthermal 
particles arranged in a power-law tail (dot-dashed). 

The high-energy tail in the ion spectrum (upper panel) 
grows with time as more and more particles are shock- 
accelerated. At late times (red curve for uj p [t = 5062), 
the slope approaches —2, which corresponds to equal en- 
ergy contributions by each decade in Lorentz factor. It 
follows that, as the upper cutoff of the spectrum increases 
linearly in time (red line in panel (b)), the energy con- 
tent of the ion nonthermal tail steadily grows, reaching 
~ 30% of the total ion energy at oj pi t = 5062. By num- 
ber, the tail is still dominated by the lower energies, and 
its fractional contribution to the ion census saturates at 
~ 5%. At late times (e.g., u) pi t = 5062), the tail does not 
resemble a simple power-law with constant slope, but it 
seems to be somewhat concave, which may be due t o th e 
modification of the shock by accelerated ions (see 93. 2|), 
as discussed for non-rela tivistic shocks in SNRs by |Am-| 
ato & Blasil (|2005l [20061). 



2200 a) 



As more ion energy is stored in the nonthermal tail, 
the peak of the ion thermal distribution shifts to lower 
energies. A similar trend is observed for electrons (lower 
panel in Fig. [7]), with the thermal bump moving to lower 
Lorentz factors as the acceleration efficiency increases]^] 
However, the tail in the electron spectrum contains less 
particles (~ 2%) and energy (~ 10%) than the ion tail, 
it is much steeper (with slope — 3.5±0.1 at uj pi t — 5062), 
and its upper energy cutoff does not grow very fast in 
time (blue line in panel (b) , multip lied by 5 for clarity) . 
Most importantly, as discussed in §3.2[ no electrons are 
seen to propagate back upstream from the shock, as op- 
posed to the large population of returning ions (compare 
the phase spaces in panels (e) and (f) of Fig. [2]). This 
suggests that there should be a substantial difference be- 
tween the acceleration process of ions and electrons, as 
we now describe. 

4.1.2. Energization Mechanism 

The orbit of a representative high-energy ion from the 
simulation of a subluminal shock (9 — 15°) is plotted 
in Fig. [8] For u p it < 1200, the particle is approach- 
ing the shock from upstream with bulk x-momentum 



-15. From uj pl t ~ 1200 to uj pi t 



2100, 



it stays at the shock gyrating around the mean magnetic 
field (see the ion x-location relative to the shock in panel 
(b)), and its Lorentz factor increases up to 7* ~ 600 
(panel (a)). Energy gain occurs primarily when the ion 
is in the upstream region, whereas its Lorentz factor does 
not significantly change while downstream (e.g., from 
Wpit ~ 1700 to uj p it ~ 1850)0 As the particle energy 
grows, its Larmor radius also proportionally increases 
(see the particle trajectory in panel (d)). The selected 
ion is finally transmitted downstream at cjpji ~ 2100, 
where it will populate the nonthermal tail seen in the 
ion spectrum of Fig. [7} In contrast, the shock- accelerated 

8 We note that, whereas the partition of energy between the ther- 
mal peak and the nonthermal tail changes with time, the average 
electron energy (that includes both components) stays relatively 
constant, at ~ 30% of the initial ion energy (blue line in Fig.ma)). 

9 This is a consequence of our choice for the simulation frame, 
which coincides with the downstream plasma frame. It follows 
that, in the downstream region, no motional electric field is present, 
and the particle energy does not appreciably change. In the shock 
frame, energy change would be seen on both sides of the shock. 
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Fig. 8. — Trajectory of a representative high-energy ion extracted 
from the simulation of a 8 = 15° subluminal shock. Panel (a): 
time evolution of the particle energy; different colors show which 
component of the electric field governs the energy change (blue for 
-E x , green for E y , red for E z ), for the portions of the trajectory 
with the largest rate of acceleration or deceleration. Panel (b): 
rr-location of the particle relative to the shock, superimposed over 
y-averaged profiles of E z stacked in time. The electric field is still 
measured in the wall frame, but shifted along x so that the shock 
appears stationary. Panel (c): time evolution of the particle 4- 
velocity, with the same color coding as in panel (a). Panel (d): 
time evolution of the particle position (x-location is relative to the 
first encounter with the shock), with the same color coding as in 
panel (a). We remind that in this case (ip = 90°) the upstream 
background electric field Eq is along — y. 



ions that are eventually reflected back upstream will con- 
tribute to the beam of returning ions seen in Fig.j2te), 
which are respo nsibl e for triggering the growth of Bell's 
instability (see £3.2 1. In turn, the waves generated by 
the instability mediate the process of ion acceleration, as 
we now describe. 

In panel (a), different colors show which component of 
the electric field governs the energy change. We see that 
energy gain is primarily associated with E y (green) and 
E z (red), which are indeed the two components associ- 
ated with Bell's circularly-polarized waves (see Fig.[3jg)). 
From the ID profiles of E z time-stacked in panel (b), 
we see that the wavelength of Bell's modes is typically 
smaller than the ion Larmor radius (compare with blue 
line in panel (b)). It follows that the interaction between 
the selected ion and the waves will be non-resonant, i.e., 
during the half Larmor period when the ion is in the 
upstream, it will encounter multiple wave fronts, result- 
ing in repeated accelerations and decelerations on a sub- 
Larmor scale. This is clearly shown in the energy evolu- 
tion of panel (a) and in the 4- velocities of panel (c) , where 
Bell's waves cause short-scale wiggles superimposed over 
the (otherwise smooth) Larmor gyration. 

The overall energy gain results from favorable encoun- 
ters (i.e., with the ion velocity locally aligned with the 
wave electric field) being more frequent than unfavor- 
able ones. In particular, significant energization may oc- 
cur when the ion is crossing the shock from upstream 
to downstream, as seen at Lo pl t ~ 1675 and u) pl t ~ 2075 
(compare panels (a) and (b)). In this case, the ion can 
remain in phase with the waves for a larger fraction of its 
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orbit, since both the ion and the waves are moving toward 
the shock. Efficient acceleration results if the wave elec- 
tric field is aligned with the ion velocity (mostly oriented 
along +Z in this phase, for a positively-charged particle 
gyrating around the mean magnetic field). Overall, the 
stochastic character of the energization process points 
towards a Fermi- like acceleration mechanism (or DSA), 
driven by non-resonant interactions of the ion with the 
self-generated Bell's waves p*] 

During the acceleration process, the ion gyrocenter 
moves predominantly along the mean magnetic field (in 
the xz plane, for tp = 90°), as shown by the time evolu- 
tion of Xi (blue) and Zi (red) in panel (d). In addition, a 
clear drift can be seen in t/j (green) towards negative val- 
ues, i.e., in the direction of the background Eo (see SS09 
for an analogous effect for positrons in pair shocks). 

As explained in p] particles reflected by the shock can 
be accelerated by the background motional electric field 
Eo while drifting along the shock surface. Here, we do 
not differentiate between SDA and SSA, but we generi- 
cally term this energization mechanism as "direct accel- 
eration" (as opposed to "diffusive acceleration") or "Eq- 
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driven acceleration". We remind that En 



Br 



is the motional field at injection, and it does not include 
any contribution from self-generated turbulence. 

In Fig. [9] we quantify the relative importance of this 
Eo-driven acceleration mechanism with respect to the 
DSA process discussed above. For each of the high- 
energy ions extracted from the simulation, we measure 
the total energy gain A-ji and the drift Axeo along the 
upstream field Eo. For each bin in A-fi, the average 
value of Axeo is shown as a black dot with error bars. 
From Axeo, we can compute the expected energy gain 
A7E0 = (<z/ mc2 ) Eq Axeo due to direct acceleration by 
the E field We sec that the resulting A7E0 (solid 
cyan line) is much lower than the actual energy gain 
(black dots), which suggests that DSA in Bell's waves 
plays a major role in the energization of ions. In fact, 
E -driven acceleration contributes at most ~ 50% of the 
total energy gain (see the dashed cyan line, which is the 
predicted A7E0 if the field Eo were twice as large). 

A separate argument needs to be made for electrons, 
whose acceleration path cannot be exactly the same as for 
ions, given the different spectrum (Fig. u\ , and above all 
the absence of returning electrons (compare phase spaces 
for ions and electrons in panels (j) and (k) of Fig. p5l re- 
spectively). In Fig. 10 we plot the orbit of a representa- 
tive high-energy electron in a 9 = 15° subluminal shock. 
The upstream trajectory of the electron gets perturbed 
(see panels (b) and (d) at uj p it ~ 2000) even before en- 
tering the shock, by the density spike formed ~ 50 c/w p j 
ahead of the shock by the pressure of returning ions. The 
electron energy (panel (a)) does not appreciably change 
before reaching the shock at u; p ;t ~ 2050. Hereafter, the 

10 We remark that the detailed nature of the waves is not instru- 
mental in the process of ion acceleration, since the Fermi mech- 
anism only requires a sufficient power in turbulent fluctuations, 
regardless of their origin. Yet, the details of wave-particle inter- 
actions (e.g., resonant versus non-resonant) may be important for 
the rate of acceleration. 

11 This formula assumes that most of the energy gain occurs in 
the upstream medium (where the background field is Eo), which 
is correct since in the simulation frame no motional electric fields 
should persist in the downstream region. 
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Fig. 9. — Relative contribution of different acceleration mecha- 
nisms, for high-energy ions (such that 7^/70 > 4 at w p jf = 2250) in 
a 8 = 15° subluminal shock. Here, Axeo is the ion displacement 
along Eo, and A-ji is the overall change in Lorentz factor. The 2D 
histogram shows the ion phase space density N(A^i, Ax^q). For 
each bin in A7i, the value of Ax^o averaged over the ion distribu- 
tion is shown as a black dot with lcr error bars. The contribution 
A7E0 expected from direct acceleration by Eo is the solid cyan 
line. The dashed cyan line is the predicted A7E0 if the field Eo 
were twice as large. 
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Fig. 10. — Trajectory of a representative high-energy electron 
extracted from the simulation of a 9 = 15° subluminal shock. Pan- 
els as in Fig. [8] but here the fluid quantity plotted in panel (b) 
is the particle number density. Also, the electron Lorentz factor 
and dimensionless momentum have been divided by the mass ratio, 
although not indicated in the plot labels. 



electron bounces several times across the shock transition 
layer (from u pi t ~ 2050 to uj pi t ~ 2175), gaining energy 
by scattering off the Bell's waves (whose wavelength im- 
prints the short-scale modulation seen in the 4- velocities 
of panel (c) for 2050 < u> pl t < 2125), similarly to what 
happens for ions. 

The final energy of the selected electron (^ e m e jm i ~ 
120) is much smaller than the downstream energy of the 
representative ion followed in Fig. [8] (7, ~ 600). This 
reflects the lower bulk energy (smaller by a factor of ~ 5, 
see Fig. [2](d)) with which electrons enter the shock, with 
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respect to ions. Due to their smaller Larmor radius, they 
are confined closer to the shock, and they do not pene- 
trate very far upstream, where the electric field of Bell's 
waves would be stronger (or, in the standard Fermi pic- 
ture, they do not sample the full velocity difference of the 
converging flows) . Also, since they are more strongly tied 
to the magnetic field, they leave the acceleration region 
earli er b y advection into the downstream (the electron in 
Fig. 10 stays at the shock for uj ps At ~ 125, as opposed 
to WpiTxt ~ 900 for the ion in Fig. [81), which justifies the 
scarcity of returning electrons^] In summary, electrons 
display a lower energy gain per acceleration cycle (rel- 
ative to ions), and a higher escape probability from the 
acceleration region. The combination of these two effects 
explains why in Fig. [7] electrons present a power-law tail 
much steeper than ions. 

As a side note, we point out that the density profiles 
stacked in panel (b) of Fig. 10 cle arly show the process 
of shock reformation discussea in §3.2| We see a density 
spike emerging ahead of the shock and then being ad- 
vected into the shock on a typical timescale of ~ 100 w"; 1 . 
Then, a new spike forms, at a slightly larger distance (in 
fact, the shock transition region widens with time), and 
the whole process starts again. 

4.2. Superluminal Shocks: cr ;t < 6 < 90° 
4.2.1. Time Evolution 
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follows the time evolution of the ion (upper 
panel) and electron (lower panel) downstream spectrum 
in a 9 — 75° superluminal shock. At all times, the ion 
spectrum is consistent with a purely Maxwellian distri- 
bution, i.e., nonthermal acceleration of ions is suppressed 
in superluminal shocks. The location of the thermal 
peak varies depending on the amount of ion energy trans- 
ferred to elec tron s before the shock, via the mechanism 
described in j 3.3 

The electron spectrum at early times (u) p it = 562, 
black curve) is also compatible with a single Maxwellian 
distribution, but at later times (blue for uj pl t = 2250, 
green for uj pi t = 3937), a separate component emerges 
at high energies. We interpret this component as the re- 
sult of extra heating of s ome electrons by the wakefield 
oscillations described in 1 3.3 rather than as nonthermal 



acceleration, for two main reasons. First, in any model of 
electron acceleration, we would expect a steady increase 
in the maximum electron Lorentz factor, at odds with 
what we observe in the blue line of panel (b). Also, we 
find that the downstream electron spectrum can be well 
fitted with two Maxwellians (panel (c)), with the hot- 
ter Maxwellian (dot-dashed line, with a temperature ~ 3 
times higher than the low-energy Maxwellian) account- 
ing for the high-energy component. A fit that employs a 
low-energy Maxwellian plus a power law with exponen- 
tial cutoff is less satisfactory. 

The appearance of the high-energy component is ac- 
companied by an increase in the average electron energy, 
at the expense of ions (panel (a) for u pl t < 3000; red for 
ions, blue for electrons). For uj pi t > 3000, the mean elec- 
tron energy decreases (panel (a)), and the high-energy 

12 Electrons are also attracted into the downstream by the cross- 
shock electric field _B X > 0, created by electron-ion charge separa- 
tion at the shock. We clearly see this effect for some electrons, but 
not for the particle in Fig 
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Fig. 11. — Time evolution of the downstream energy spectrum in 
a = 75° superluminal shock, for ions (upper panel) and electrons 
(lower panel): oj p ;t = 562 (black), oj p ;t = 2250 (blue), oj p ;t = 3937 
(green), cj p jt = 5625 (yellow), and w p jf = 7312 (red). Subpanels 
as in the caption of Fig. [7] but here the electron spectrum in panel 
(c) is fitted with two Maxwellians. 



component recedes (wpjf = 5625, yellow curve), until it fi- 
nally disappears (u) p it = 7312, red curve). The similarity 
between the time evolution of the average electron energy 
in panel (a) and the maximum electron Lorentz factor in 
panel (b) (a good proxy for the importance of the high- 
energy component) suggests that a common mechanism 
should regulate both the low-energy "thermal" part and 
the high-energy tail, as we now describe. 

At the beginning, electrons enter the shock with x- 
momentum = — 7o/?o — — 15 and a small thermal spread. 
The synchrotron maser instability at the shock front pro- 
duces powerful electromagnetic waves (Poynting flux in 
excess of cEfjAit) that illuminate the upstream flow, 
causing efficient t rans fer of bulk momentum from ions 
to electrons (see j ]3.3[ ). As electrons with larger bulk 
energy enter the shock, the synchrotron maser insta- 
bility is more efficient in producing strong precursor 
waves, and the transfer of energy from ions to electrons 
in the upstream becomes even more significant. This 
self-reinforcing cycle repeats until electrons reach energy 
equipartition with ions (panel (a), at uj pl t ~ 2500). As 
the strength of the electromagnetic precursor increases, 
wakefield oscillatio ns of larger am plitude are generated 
in the upstream ( |Lyubarsky| |2006j ) . As anticipated in 
|3.3| their dissipation then results in a more pronounced 
high-energy component in the electron spectrum. This 
explains why the maximum electron Lorentz factor in- 
creases for w p i£ < 3000 (blue line in panel (b)). 

Heated by the upstream wakefield oscillations, elec- 
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trons now enter the shock much hotter than before, i.e., 
their comoving momentum spread is larger. In such con- 
ditions, th e growth of the synchrotr on maser instability is 
quenched ( Hoshino fc Arons||1991 ), and the electromag- 
netic precursor becomes less powerful. This reduces the 
transfer of bulk energy from ions to electrons in the up- 
stream, causing the decrease in avera ge e lectron energy 
observed at Lu pi t > 3000 in panel (a)Fj Also, the up- 
stream wakefield oscillations now become weaker, so that 
fewer electrons are injected into the high-energy spectral 
component, and the maximum electron Lorentz factor 
decreases (panel (b), at oj pl t > 3000). 

In our simulation, we find that the electron spectrum 
approaches a steady state for uj p [t > 7000 (we followed 
the shock evolution up to w p ;i ~ 10, 000). However, since 
the incoming electrons are now colder, we cannot exclude 
that at later times the whole process may repeat, with 
efficient transfer of energy to electrons accompanied by 
the emergence of a new high-energy spectral component. 
Instead, we do not expect the ion spectrum to deviate 
from a Maxwellian distribution, since ions are not signif- 
icantly affected by the upstream wakefield oscillations, 
due to their larger inertia. 

Finally, we comment on the difference between our 
2D sim ulations and the ID results presented by |Hoshino| 
(2008), that reported nonthermal acceleration of elec- 
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trons in perpendicular magnetized electron- ion shocks. 
First, we find that ID simulations over-estimate the 
power of the electromagnetic precursor, since they can- 
not resolve variations along the shock surface (Fig. [6^c) 
shows that the strength of the electromagnetic precursor 
is not uniform in y). As a result, ID runs tend to over- 
estimate the efficiency of ion-to-electron energy transfer, 
and the importance of the high-energy component in the 
electron spectrum. For nii/m e = 16, computational do- 
mains with more than 64 transverse cells are requ i red to 
obtain consistent 2D results^ Second, Hoshino (20081 
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followed the shock evolution only up to oj p it ~ 850, a 
timespan too short to appreciate the saturation and de- 
crease in the maximum electron Lorentz factor shown in 
panel (b). As mentioned above, this piece of evidence is 
essential for our interpretation of the high-energy spec- 
tral component as a result of wakefield heating, rather 
than acceleration. 

4.2.2. Energization Mechanism 

The trajectory of a representative high-energy elec- 
tron in a superluminal shock (9 = 75°) is presented in 
Panel (a) shows that all of the energy gain, up 
70, occurs in the upstream region before 
the electron encounters the shock at uj p it ~ 1050 (see the 
electron x- location relative to the shock in panel (b)). 
After that, the electron is advected downstream, with no 
appreciable change in energy. 

13 This effect has been verified with controlled numerical exper- 
iments, in which we initialized the incoming particles with a larger 
thermal spread than the usual A7 = 10 — 4 . For A7 > 10 — 1 , both 
the power of the electromagnetic precursor and the upstream ion- 
to-electron energy transfer are gradually suppressed, as well as the 
high-energy component in the electron spectrum. 

14 In fact, to follow the spectral evolution shown in Fig. |ll| we 
have used a simulation box with 64 transverse cells, which is small 
enough to suppress the growth of the numerical instability men- 
tioned in ij2] yet large enough to preserve the relevant 2D properties 
of the shock. 




Fig. 12. — Trajectory of a representative high-energy electron 
extracted from the simulation of a 9 = 75° superluminal shock. 
Panels as in Fig.[8] but here the fluid quantity plotted in panel (b) is 
-E x . Also, the electron Lorentz factor and dimensionless momentum 
have been divided by the mass ratio, although not indicated in the 
plot labels. We remind that in this case (ip = 0°), the upstream 
background electric field Eq is along +z. 



For u> p it < 850, the electron is first boosted and then 
deboosted (see energy, panel (a); or x- momentum, blue 
line in panel (c)), in response to the wakefield wiggle 
seen at ~ 350c/u;pi ahead of the shock (E x in panel (b)). 



So far, the trajectory of the electron in Fig. [12] presents 
no qualitative difference from the bulk of incoming elec- 
trons. But at ujpit ~ 850, presumably due to a strong 
upstream-oriented kick imparted by the wakefields, the 
energy of the selected electron becomes so small that it 
can eas ily decouple from the bulk flow. As anticipated 
in S3.3 this is the "injection" step for all the electrons 
that will end up in the high-energy spectral component 
seen in Fig. |TT) 

After decoupling from the bulk flow, the electron gy- 
rocenter starts moving with drift velocity E x B /-Bg, 
which for the field configuration employed here (cp = 0°) 
corresponds to a straight path in the xy plane (see y e in 
panel (d) for uj p it > 850). At the same time, the electron 
gets pulled by the upstream motional electric field Eq, 
and it performs a series of Larmor cycles (see energy in 
panel (a) and 4- velocity in panel (c), for uj pi t > 850). 
When the electron velocity is aligned with — Eo the par- 
ticle energy grows, but it decreases by the same amount 
when the two are oppositely directed (see red segments in 
panel (a)). Overall, no net energy gain is expected from 
Eo during a complete Larmor cycle. In fact, as seen from 
the upstream frame, the electron motion is just a Larmor 
gyration around the background magnetic field, with en- 
ergy determined by the magnitude of the kick imparted 
by the wakefields. 

The final electron energy is determined by the combi- 
nation of two factors: (1) the amplitude of the energy 
oscillations seen in panel (a); and (n) the gyro-phase 
with which the particle enters the shock. Regarding (i), 
secular variations (i.e., averaged over gyro-phase) in the 
electron energy are determined by the kicks imparted by 
the wakefields. In particular, the electron will be most 



14 



sensitive to such kicks when its energy is smaller. At this 
point, as seen from the upstream frame, the electron is 
moving away from the shock. If E K < 0, the electron 
will receive a kick in the same direction of its momen- 
tum, and its upstream-frame energy will increase. As 
seen from the downstream frame, this corresponds to a 
larger amplitude of the resulting energy oscillation (e.g., 
at x ~ 960c/a;pi in panel (a)). The opposite happens if 
the phase of wakefield oscillations is such that E x > 0, so 
that the following energy oscillation will have a smaller 
amplitude (e.g., at x ~ 900c/w p ; in panel (a)). As a con- 
firmation of the importance of wakefields, we see that the 
amplitude of the energy oscillations in panel (a) increases 
dramatically for uj pi t > 925, when the selected electron 
enters a region with strong wakefields (E x in panel (b), 
following the particle track for u! p [t > 925). 

On the other hand, depending on the gyro-phase 
with which the particle encounters the shock (point (ii) 
above) , the electron will end up in a different energy bin 
of the downstream high-energy Maxwellian (dot-dashed 
line in Fig. [FIJc)). 

Finally, we remark that, despite the significant wake- 
field heating experienced in the upstream, which could 
possibly serve as pre-injection for a Fermi-like process, 
the electron in Fig. [12] does not bounce back from the 
shock into the upstream, due to insufficient turbulence on 
the downstream side. In fact, in our simulations of super- 
luminal shocks, we do not find any evidence for return- 
ing electrons (see the electron phase space in Fig. [5jf ) ) , 
confirming that Fermi acceleration is suppressed for su- 
perluminal configurations. 

5. DEPENDENCE ON UPSTREAM CONDITIONS 

Armed with a better understanding of the energiza- 
tion mechanisms that operate in relativistic magnetized 
electron-ion shocks, we now explore how the efficiency 
of ion-to-electron energy transfer and of particle accel- 
eration depends on the conditions of the upstream flow. 
First, keeping the bulk Lorentz factor (70 = 15) and 
magnetization (<r = 0.1) fixed, as we did so far, we in- 
vestigate the transition between subluminal and super- 
luminal shocks by exploring the full ra nge of magnetic 
obliquities, from 9 = 0° to 9 = 90° (in [O). Then, for 
the two representative obliquities discussed above (one 
subluminal, 9 = 15°; one superluminal, 9 = 75°), we ex- 
tend our analysis to different values of the magnetization 
(in 1 5.2 1 and bulk Lorentz factor (in j |5.3| ). 

The comparisons in this section are performed at 
ajpii = 2250. Yet, as shown in Figs. [7] and 11 the down- 
stream particle spectrum is not in steady state. In sub- 
luminal shocks (Fig. m), the nonthermal tail grows with 
time, so that the acceleration efficiency we compute at 
ujpit — 2250 should be t aken as a lower limit. In super- 
luminal shocks (Fig. 11), the electron high-energy com- 



ponent is most pronounced between u> p it ~ 2000 and 
Upit ~ 3000, and then it disappears, althoug h a cy clical 



behavior is likely, yet undemonstrated (see S 4.2.1 1. At 
uj pl t — 2250, we are capturing this high-energy compo- 
nent when it is close to its maximum, and more likely to 
result in an observational signature. 

In the following, we fit the downstream particle spec- 
trum with a three-dimensional Maxwellian plus a power- 
law tail with an exponential cutoff. This is not strictly 
appropriate for superluminal magnetized shocks, where 



the high-energy component of the electron spectrum is 
better fitt ed wit h a second (hotter) Maxwellian, as dis- 
cussed in §4.2.1| Although this component actually re- 
sults from heating (rather than acceleration) , in this sec- 
tion we treat it as a nonthermal tail, for easier compari- 
son with observations F^l 

The form chosen for our fitting function allows to split 
the particle population into a thermal and a nonthermal 
component, and to follow how the relative importance 
of the two sub-popu lations changes with the upstream 
parameters. In Figs. 13p7 panel (a) presents the mean 
downstream particle energy (red for ions, blue for elec- 
trons) in units of the bulk energy of injected ions, in- 
cluding both thermal and nonthermal particles. Panels 

(b) -(d) provide a better characterization of the nonther- 
mal component: panel (b) shows the slope p of the best- 
fitting power law, such that dN Ph (j) / d~/ cx 7 _p ; panels 

(c) and (d) present respectively the fraction of particles 
and energy contained in the nonthermal tail (red for ions 
and blue for electrons, relative to the total number or en- 
ergy of that species)^] In subpanels (b)-(d), dotted lines 
(instead of the usual solid lines) mark the cases in which a 
nonthermal power-law tail is not a physically-motivated 
model for the high-energy component (i.e., for electrons 
in superluminal magnetized shocks, as discussed above). 

5.1. Dependence on Magnetic Obliquity 9 

Fig. [13] presents the downstream particle spectrum at 
ujpit — 2250 (upper panel for ions, lower panel for elec- 
trons), for a sample of magnetic obliquities covering the 
full range from 9 = 0° to 90°. In this section, we de- 
scribe how the efficiency of particle acceleration and elec- 
tron heating depends on magnetic inclination, emphasiz- 
ing the importance of the critical boundary cl .[ t ~ 34° 
between subluminal and superluminal configurations. 

The ion spectrum in subluminal shocks (blue curve for 
9 = 0°, green for 9 = 15°, red for 9 = 30°) shows a promi- 
nent nonthermal tail at high energies, together with a 
thermal bump at low energies. Instead, for superluminal 
shocks (black for 9 — 45°, yellow for 9 = 75°, purple for 
9 = 90°), the ion spectrum does not appreciably devi- 
ate from a Maxwellian distribution. As shown in panels 

(c) and (d) , the efficiency of ion acceleration (red lines) 
abruptly drops between 9 = 30° and 9 = 35°, as the 
magnetic obliquity passes the critical value 9 clit ~ 34°. 

For subluminal configurations, the ion nonthermal tail 
shows a non-monotonic behavior with angle. Ion acceler- 
ation is most efficient for 9 — 15°, with a fraction ~ 4% 
of particles and ~ 22% of energy contained in the tail 
at ujpit = 2250 (red dots for 9 = 15° in panels (c) and 

(d) , respectively). As discussed in £4.1.1 at uj p \t = 5062 
the ion acceleration efficiency reaches as much as ~ 5% 
by number and ~ 30% by energy (red crosses in panels 
(c) and (d)). By comparison, at uj pi t = 2250 the ion 
nonthermal tail for 9 = 0° only accounts for ~ 2% of 
particles and ~ 10% of ion energy, whereas the accelera- 

15 Observationally, a suprathermal component in the radiation 
spectrum will most likely be ascribed to nonthermal acceleration 
of particles, and fitted as a power-law. 

16 For instance, to obtain the fraction of energy in nonthermal 
electrons relative to the bulk energy of injected ions (which is often 
required in phenomenological models of GRBs and AGN jets), the 
value read from the blue line in panel (d) should be multiplied by 
the value of the blue line in panel (a). 
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tion efficiency for 9 = 30° (~ 4% by number, ~ 20% by 
energy) is comparable to 9 — 15°, but the tail extends 
to smaller Lorentz factors. Also, the ion tail is flatter for 
9 = 15°, with a slope p ~ 2.2, compared to p ~ 3.5 for 
9 = 0° and p ~ 2.5 for 9 = 30° (red line in panel (b)). 

At earlier times (e.g., uj pi t = 1125, not shown here), 
the trend was different, with larger angles (yet still sub- 
luminal) showing flatter ion tails that extended to higher 
Lorentz factors and contained more particles and energy. 
In fact, this was the hierarchy that SS09 observed in 
electron-positron shocks (see Fig. 11 in SS09, same col- 
ors as here). In SS09, such trend was attributed to a 
change in the dominant acceleration mechanism, with 
DSA prevailing in quasi-parallel shocks and SDA most 
efficient for 15° < 9 < 9 C Iit . In electron-ion shocks, as 
we have shown in §4.1.2t magnetic turbulence created 
by the returning ions (in the form of Bell's waves) can 
dramatically affect the acceleration process, much more 
than in electron-positron shocks, where the amplitude of 
self-generated upstream turbulence was much lower (see 
SS09). T he h ierarchy seen in the subluminal ion spec- 
tra of Fig. [13] will then reflect the influence of magnetic 
inclination on the strength of Bell's waves. 

Since the gyrocenters of returning particles tend to 
move along the upstream field, for higher (yet sublu- 
minal) obliquities the returning ions will be confined 
closer to the shock. This has two opposite consequences. 
First, since the returning ions are responsible for trig- 
gering Bell's instability, at larger obliquities its growth 
is more likely to be suppressed by advection into the 
shock. In fact, at u p it = 2250, Bell's modes are already 
fully developed for 9 — 15°, whereas for 9 = 30° they 
are just starting to grow, and the acceleration process 
is still governed by SDA, as in electron-positron flows 
(SS09). On the other hand, Bell's instability saturates 
when the self-generated magneti c energy approaches the 



energ y density of returning ions ( Riquelme & Spitkovsky 
2009 1 , which is larger for higher obliquities since the re 



turning particles accumulate in a smaller region ahead of 
the shock, as discussed above. It follows that Bell's waves 
are stronger for 9 = 15° than for 6 — 0°. Overall, the 
best compromise at cj p ;t = 2250 is achieved for 9 = 15°, 
where the waves are both powerful and rapidly-growing 
(compared to advection). 

At later times, the situation may change. The growth 
of Bell's waves for 9 = 30° will not be limited by ad- 
vection any longer, and the turbulence will saturate at 
higher amplitudes than for 9 = 15°, giving efficient 
ion acceleration via DSA. Therefore, the trend in ac- 
celeration efficiency observed for subluminal electron- 
positron shocks by SS09 may eventually be restored for 
ions in electron-ion shocks as well. However, here it 
will be driven by the increase with angle in the self- 
generated turbulence that mediates DSA, whereas in 
electron-positron flows it follows the amplitude of the 
background motio nal e lectric field responsible for SDA. 

As discussed in |4.1| acceleration of electrons in sub- 
luminal shocks is depressed with respect to ions. For 
nearly-parallel shocks (blue curve for 9 = 0°, green for 
9 = 15°), a minor fraction (~ 1%) of the incoming elec- 
trons are accelerated during their passage through the 
shock, as explained in { 4.1.2| In the downstream, they 
populate a steep (p ~ 5) power-law tail which only con- 
tains ~ 3% of the total electron energy. The tail may get 
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Fig. 13. — Downstream particle spectrum at u> p it = 2250 (up- 
per panel for ions, lower panel for electrons) for different magnetic 
obliquities, from 8 = 0° (parallel shock) to 6 = 90° (perpendicular 
shock). In subpanels, subluminal and supcrluminal geometries are 
separated by a vertical dotted black line, corresponding to the crit- 
ical angle 6 cr i t ~ 34°. Subpanels: (a) downstream ion (red) and 
electron (blue) average energy, in units of the bulk energy of in- 
jected ions, with the horizontal dashed blue line showing the value 
expected for electrons in the absence of any ion-to-electron energy 
transfer (= m e /mi ~ 0.06); (b) power-law slope of the nonthermal 
tail (red for ions, blue for electrons); (c) fraction of ions (red) or 
electrons (blue) in the nonthermal tail; (d) fraction of energy in 
the ion (red) or electron (blue) nonthermal tail, with respect to 
the total kinetic energy of that species. In subpanels (b)-(d), dot- 
ted blue lines mark the cases in which a nonthermal power-law tail 
is not a physically-motivated model for the electron high-energy 
component (see text). Crosses for 9 = 15° (red for ions, blue for 
electrons) correspond to values measured at ui p it = 5062 (instead 
of ojpit = 2250), to show how the nonthermal tail changes in time. 



somewhat flatter with time, thus accounting for more 
particles and energy (see blue crosses in panels (b)-(d), 
for 9 = 15°), but the acceleration efficiency for electrons 
remains smaller than for ions by at least a factor of three. 
Close to the critical obliquity 9 CI i t — 34° (see 9 = 30°, 
red curve), a more pronounced tail appears in the elec- 
tron spectrum, populated by particles accelerated at the 
shock via the SDA mechanism, in analogy to what hap- 
pens in electron-positron flows (SS09). However, despite 
extending to high energies (up to 7 e m e /mi ~ 400 at 
u p it = 2250), the electron tail for 9 = 30° only contains 
~ 1% of particles and ~ 2% of electron energy. 

Instead, a prominent high-energy tail is observed in 
the electron spectrum of superluminal shocks, beyond 
the thermal bump. As the obliquity angle increases from 
9 = 75° (yellow curve) to 9 = 90° (purple), the tail 
stretches to higher Lorentz factors (maximum 7 e m e /roi 
from ~ 95 to ~ 150), and it contributes a larger fraction 
of electrons (from ~ 2% to ~ 13%) and electron energy 
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(from ~ 6% to as much as ~ 35%). If this high-energy 
component is fitted with a power law (but a fit with a 
seco nd Ma xwellian is actually more physically motivated, 
see f 4.2.1 ), we find that it gets flatter for higher obliqui- 
ties (Ifom p ~ 2.3 fo£0 = 75° top - 1.5 for 9 = 90°). 

As explained in §4.2. 2| this high-energy component re- 
sults from extra heating of some electron s by the up- 
stream wakefield oscillations described in 93.31 that are 
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ultimately generated by the radiative push of the pre- 
cursor wave on the incoming flow. Since the amplitude 
of wakefield oscillati ons increases wit h the strength of 
the precursor wave ( jLyubarsky 2006 ) , which is larger 
for higher obliquities (roughly oc Efi oc sin 2 0), the elec- 
tron high-energy component will be more pron ounced for 
quasi-perpendicular shocks, as seen in Fig. |13[ Being due 
to heating, rather than acceleration, we do not expect the 
electron tail in superluminal shocks to extend with time 
to higher Lorentz factors. Rather, it may recede to lower 
energies, and possibly be hidde n by the low-energy ther- 
mal bump, as shown in { 4.2.1 for 9 — 75°. 

In superluminal shocks, the electromagnetic precursor 
wave also causes efficient transfe r of energy from ions to 
electrons ahead of the shock (see { 3.3 ). As a result, in the 
downstream region the el ectr on and ion average energies 
are comparable (see Fig. [l3j[a), for 35° < 6 < 75°), or 
electrons may be even hotter than ions (for 75° < 9 < 
90°). In subluminal shocks (0° < 9 < 30°), electrons are 
boosted and heated ahead of the shock, as a result of the 
perturbation induced on t he incoming flow by the cloud 
of returning ions (see j j3.2| . With respect to superluminal 
angles, here the downstream electron energy is a smaller 
fraction (~ 20% — 30%) of the initial ion bulk energy, but 
still much in excess of m e /m,i ~ 6% (for rrii/m e = 16), 
the value we would expect in the absence of any ion-to- 
electron energy transfer. 

In Appendix [3] we show that the values quoted here 
for the efficiency of electron heating and of particle ac- 
celeration are essentially unchanged for larger (and more 
realistic) mass ratios. 

5.2. Dependence on Magnetization a 

In this section, we explore the dependence of shock 
thermalization and acceleration upon the magnetization 
of the upstream flow. We keep the magnetic obliquity 
fixed and vary the magnetization from a = 10~ 5 , a vir- 
tually unmagnetized shock, up to a = 1.0, where the 
upstream field is in equipartition with the kinetic energy 
of injected part icles. We discuss separately the cases of 
sublum inal (in \ 5.2.1 for 6 — 15°) and superluminal (in 
Sj5~Z2) for 8 = 7V7shocks. 



5.2.1. Subluminal Shocks: 0° < 6 < cr i t 

Fig. [14] presents the downstream particle spectrum of 
a 9 = 15° subluminal shock, for different magnetiza- 
tions. The ion spectrum (upper panel) shows in all cases 
a pronounced nonthermal tail, whereas a significant high- 
energy component appears for electrons (lower panel) 
only at a < 10~ 2 . This corresponds to a change in the 
mechanism that mediates the shock. 

For a < 10~ 2 (blue curve for a = 10 -2 , red for 
a = 10~ 3 ), the shock is mediated by the filamentation 
(Weibel) instabili t y, as i n strictly unmagnetized flows 
(Spitkovsky 2005 2008a I. The free energy for the in- 
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Fig. 14. — Downstream particle spectrum at u> p it = 2250 (upper 
panel for ions, lower panel for electrons) for different magnetiza- 
tions. We fix the m agne tic obliquity 9 = 15° (subluminal shock). 
Subpanels as in 
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stability comes from the counter-streaming between the 
incoming plasma and the shock-accelerated particles that 
propagate ahead of the shock. In turn, magnetic turbu- 
lence generated by the instability mediates nonthermal 
acceleration of p articles at the shock front, via a Fermi - 
like mechanism ( |Spitkovsky||2008b| partins et al.||2009[) . 

In fact, the ion spectrum m such shocks shows a promi- 
nent nonthermal tail. As the magnetization decreases be- 
low a = 10~ 2 , the tail becomes flatter (from p — 3.5 at 

2 ^ O r, „f „ _ 1 D-3 raA lino in TTirr Fj"^ ' \X\ 

5% 
to 



a = 10~ 2 to p — 2.5 at o = 10~ 3 , red line in Fig. 
and it contains a larger fraction of ions (from 
to — 7%, Fig. fliTc)) and ion energy (from — 17% 
- 24%, 
o- < 10 



Fig. 14 d)). Similarly, the electron spectrum for 



2 presents a significant high-energy tail, compa- 
rable to the ion tail in slope (compare red and blue lines 
in panel (b), for a < 10~ 2 ) and maximum Lorentz fac- 
tor. The acceleration efficiency for electrons is not much 
smaller than for ions, at most a factor of two (compare 
red and blue lines in panels (c) and (d), for a < 10~ 2 )FH 
The similarity between the acceleration properties of 
electrons and ions in weakly magnetized flows (a < 10~ 2 ) 
comes from the fact that the two species enter the shock 
with comparable bulk energy (compare red and blue 
lines in panel (a), for a < 10~ 2 ), due to efficient trans- 
fer of energy from ions to electrons ahead of the shock 
(Spitkovsky et al., in prep). The Fermi mechanism, 

17 For both electrons and ions, the downstream spectrum at 
even lower magnetizations (we explored down to a = 10 — 7 ) is 
identical to the case a = 10 -3 shown here, which is therefore a good 
representation of weakly magnetized (and unmagnetized) shocks. 
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which only depends on the particle Larmor radius (i.e., 
on particle energy, rather than mass), will then proceed 
in a similar way for ions and electrons. 

The situation changes for a > 10~ 2 . Here, the mag- 
netic field is strong enough so that the upstream plasma 
cannot filament, and the growth of Weibel instability is 
suppressed. In these conditions, the current of shock- 
accelerated ions moving ahead of the shock can trigger 
Bell's instability, as discussed in £3.2 For 10~ 2 < a < 
0.3, the waves generated by Bell's mechanism are respon- 
sible for mediating the shock and driving the acceleration 
of ions. They are most powerful for a ~ 0.1 (black curve), 
which then corresponds to a peak in the ion acceleration 
efficiency (~ 4% by number, ~ 22% by energy) and a 
flatter nonthermal tail (p ~ 2.2). 

In contrast, acceleration of electrons is depressed in 
highly ma gneti zed flows (black curve for a — 0.1), as dis- 
cussed in §4.1| As the magnetization increases, return- 
ing particles are confined closer to the shock (in units 
of c/wpi, but same distance in Larmor radii), so their 
perturbation to the incoming flow extends over a smaller 
length. It follows that the transfer of energy from ions to 
electrons in the upstream has less time to proceed, and 
it is therefore less efficient (see panel (a), from a = 10~ 3 
to (7 = 0.1). For a — 0.1, electrons e nter the shock with 
~ 15% of the initial ion energy (see |3.2[ ), so their Lar- 
mor radius is small compared to ions. Then, only a small 
fraction of electrons (~ 1%, panel (c)) are accelerated at 
the shock before being advected downstream, where they 
populate a steep nonthermal tail {p ~ 4.5, panel (b)). 

At even higher magnetizations (a > 0.3, yellow curve 
for a = 0.3), the instability that mediates the shock 
changes again, and the dominant mode is now generated 
via gyro- frequency resonance with the returning ions (see 
Appendix [C). However, the overall shape of the electron 
and ion spectrum does not significantly change (compare 
black, for a — 0.1, and yellow, for a — 0.3), showing 
that, despite the different polarization of the dominant 
mode (now resonant, as opposed to Bell's non-resonant 
waves) , the physics of the acceleration is not significantly 
altered. However, the level of magnetic turbulence gener- 
ated by the instability becomes smaller for a > 0.1, which 
explains the observed decrease in acceleration efficiency, 
especially for ions (Fig. [Ti^c) and (d), for a > 0.1). 

The results presented m Fig. 14 refer to uj^t = 2250. 
At later times, the nonthermal tail will extend to higher 
Lorentz factors in all cases, yet at a different rate for dif- 
ferent magnetizations. Since ion acceleration is basically 
regulated by the degree of magnetic turbulence, which is 
the largest for a ~ 0.1, the fastest evolution in the ion 
spectral tail will be observed for a ~ 0.1 (black curve), 
which in fact extends to the highest energy. Instead, the 
electron tail evolves more rapidly for moderate magne- 
tizations (a = 10~ 2 in blue), since electron acceleration 
in high-er flo ws i s impeded by advection into the down- 
stream (see §4TJ). In g ener al, the hierarchy in maximum 
energy observed in Fig. [14] as a function of a corresponds 
to an analogous trend in the rates of acceleration. 

Finally, we remark that the results reported here for 
9 = 15° are representative of the general behavior of 
low-obliquity shocks. 

5.2.2. Superluminal Shocks: f? C rit < < 90° 
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Fig. 15. — Downstream particle spectrum at u> p it = 2250 (upper 
panel for ions, lower panel for electrons) for different magnetiza- 
tions. We fix the ma gne tic obliquity 8 = 75° (superluminal shock). 
Subpanels as in 
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The downstream particle spectrum at different magne- 
tizatio ns fo r a superluminal shock with 9 = 75° is shown 
in Fig. 15 The ion spectrum (upper panel) is always 
consistent with a purely Maxwellian distribution, with 
no evidence for nonthermal particles, except for the low- 
est magnetization (a = 10 -5 , green curve). The non- 
Maxwellian shape of the ion spectrum for a = 10~ 3 (red 
curve) results from incomplete ion thermalization at the 
shock front, and it will eventually relax to a Maxwellian 
distribution further downstream from the shock. In con- 
trast, a prominent high-energy component appears in the 
electron spectrum (lower panel) for both very low (green 
for g = 10 -5 ) and very high (black for a = 0.1, yellow for 
<7 = 0.3) magnetizations, but not in between. However, 
the mechanism that regulates the electron high-energy 
tail is markedly different between low and high magneti- 
zations, as we now describe. 

In high-cr flows (a > 0.1; black curve for a = 0.1, 
yellow for a = 0.3), the electromagnetic precursor wave 
generated by the synchrotron maser instability ensures 
efficient transfer of energy fr om i ons to electrons ahead 

For a > 0.1, the pre- 



of the shock, as discussed in \ 



3.3 



cursor is strong enough to drive electrons up to energy 
equipartition with ions (compare red and blu e lines in 



panel (a) for a > 0.1). As described in ^4.2 powerful 



wakefield oscillations excited by the precursor wave in 
the upstream region dissipate their energy by producing 
a pronounced high-energy tail in the electron spectrum 
(black for a = 0.1, yellow for a = 0.3). The tail contains 
~ 2% of electrons and ~ 8% of electron energy (blue 
line in panels (c) and (d), respectively). Again, we re- 
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mind that in superluminal magnetized shocks (cr > 0.1), 
the electron high-energy tail is due to wakefield heating 
(rather than acceleration), and it d oes not extend in time 
to higher Lorentz factors (see Fig. 11 for a = 0.1). 

10~ 3 < a < 0.1; 
= 10~ 3 ), the syn- 
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At intermediate magnetizations ( 
blue curve for a = 10~ 2 , red for a 
chrotron maser instability is gradually suppressed, since 
it is harder for the electron motion at the shock to stay 
coherent on the longer gyration tim e. It follows that th e 
electromagnetic precursor is weaker ( |Gallant et al.|1992[ ), 
and ion-to-electron energy transfer in the upstream is 
suppressed. As a result, the electron spectrum shifts to 
lower energies, whereas ions get hotter (see panel (a), 
from a = 0.1 down to a = 10~ 3 ). A minor high-energy 
component, still due to wakefield heating, survives in 
the electron spectrum at a = 10~ 2 (blue curve), but 
it disappears for a = 10~ 3 (red curve). Despite the in- 
efficiency of electron heating upstream of the shock, for 
a = 10~ 3 the downstream electron population still con- 
tains as much as ~ 20% of the initial ion kinetic energy 
(blue line in panel (a)). This is considerably in excess of 
m e /mi ~ 6% (dashed blue line in panel (a)), the value 
we would expect in the absence of any ion-to-electron en- 
ergy transfer^] In this regime, electrons gain energy not 
in the upstream region, but directly at the shock front 
by the cross-shock potential in the ion foot. 

At lower magnetizations (10~ 5 < a < 10~ 3 ; purple 
for a — 10~ 4 ), the ion foot, whose scale is set by the 
Larmor radius of incoming ions in the shock-compressed 
field, becomes wider, so that electron heating by the 
cross-shock potential is more efficient (panel (a), for 
a < 10~ 3 ). Finally, for a < 10~ 5 (we have tested down 
to (j = 10~ 7 ), the ion foot is so wide that the counter- 
streaming between incoming and reflected ions can trig- 
ger the Weibel instability, before the upstream flow gets 
advecte d into the shock. As observed for unmagnetized 
shocks ( Spitkovsky 2008a), this is accompanied by strong 
electron heating m the Weibel filaments and electron and 
ion acceleration via a Fermi-like mechanism. This ex- 
plains the nonthermal tail seen in the electron and ion 
spectra of low-cr flows (green curve for a = 10~ 5 ). For 
ions, the tail is relatively flat (p ~ 2.7, panel (b)) and it 
contains ~ 3% of ions and ~ 13% of ion energy (panels 
(c) and (d), respectively). The acceleration efficiency for 
electrons is smaller by roughly a factor of two (compare 
red and blue lines in panels (c) and (d), for a = 10~ 5 ). 
Unlike in high- a shocks, here the ion and electron non- 
thermal tails extend with time to higher energies. 

For a — 10~ 5 , the ion and electron spectra in Fig. 15 
approach the results obtained for 9 = 15° at low magne- 
tizations (e.g., a — 10~ 3 , red curve in Fig. 14 1. Indepen- 
dently of magnetic obliquity, the s hock now behaves as 



the unmagnetized case discussed by Spitkovsky] ([2008a] 



We find that the transition to the unmagnetized regime 
happens at lower magnetizations for higher obliquities 
(cr - 10~ 2 for 6 = 15°, cr ~ 10~ 4 for 9 = 75°), when the 
self-generated Weibel fields are comparable to the trans- 
verse component of the background field (cx ^Ja sin 9) . 
Neglecting the obliquity dependence, a ~ 10~ 3 may be 
taken as a representative value for the transition point. 

18 This result has been confirmed with larger mass ratios, up to 
m,i/m e = 1000. 




it = 2250 (upper 

panel for ions, lower panel for electrons) for different upstream bulk 
Lorentz factors. We fix the magnetic obliquity 9 = 15° (subluminal 
shock) and the magnetization cr = 0.1. Spectra are shifted along 
the x-axis by 15/7o to facilitate comparison with the reference case 
70 = 15. 



We remark that the results presented here for 9 = 75° 
are a good representation of the overall behavior of su- 
perluminal shocks, up to the perpendicular case 9 = 90° . 

5.3. Dependence on Bulk Lorentz Factor 70 

In this section, we keep the magnetization fixed at a — 
0.1 and explore the dependence of our results upon the 
upstream bulk Lorentz factor, from 70 = 3 to 70 = 50. 
We analy ze separately the case of sublu minal (9 = 15°, 
in f 5.3.1 1 and superluminal (9 = 75°, in ( 5.3.2) shocks. 



5.3.1. Subluminal Shocks: 0° < 9 < # cr it 

Fig. [TBI shows the downstream particle spectrum of a 
9 = 15° subluminal shock, for different upstream bulk 
Lorentz factors. Spectra are shifted along the x-axis by 
I5/70, to facilitate comparison with the reference case 
70 = 15 discussed in the previous sections. Once normal- 
ized in this way, the spectra for 70 > 15 overlap almost 
perfectly (black for 70 = 15, blue for 70 = 30, red for 
70 = 50), suggesting that the physics that regulates the 
acceleration of particles (both ions and electrons) is in- 
sensitive to the upstream bulk Lorentz factor, for highly 
rclativistic flows. In fact, in our simulations we find that 
f° r 7o ^ 15 the waves generated by Bell's instability 
have the same power and wavelength, once their power 
is normalized to the upstream bulk kinetic flux and the 
wavelength is measured in units of the upstream ion skin 
depth (or Larmor radius). The scaling seen in Fig. 16 
for 70 > 15 then follows from the fact that Bell's modes 
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Fig. 17. — Downstream particle spectrum at u v \t = 2250 (up- 
per panel for ions, lower panel for electrons) for different upstream 
bulk Lorentz factors. We fix the magnetic obliquity 9 = 75° (su- 
perluminal shock) and the magnetization a = 0.1. Spectra are 
shifted along the x-axis by 15/7o to facilitate comparison with the 
reference case 70 = 15. 



are of paramount impor tance for particle acceleration in 
subluminal shocks (see [ 4.1.2 1. 

Spectra for lower Lorentz factors (green for 70 = 5, 
yellow for 70 = 3) still share the same overall properties 
of 70 = 15, albeit with a tendency for lower high-energy 
cutoffs and cooler electron distributions. Overall, the 
properties of 70 > 3 subluminal shocks seem to be well 
captured by our representative case 70 = 15. 

In mildly relativistic magnetized flows (70 < 2, for 
a = 0.1), the shock is mediated by a different process, 
such that the dominant mode is now generated by gyro- 
frequency resonance with the returning ions, rather than 
by Bell's non-resonant instability (see Appendix[C|. This 
may change the efficiency of particle acceleration. Such 
trans-relativistic regime is not the subject of this work 
and will be addressed elsewhere. 



5.3.2. Superluminal Shocks: r? cr it < < 90° 

As shown in Fig. [T7J the spectrum of 9 = 75° super- 
luminal shocks presents a much more complicated de- 
pendence on the upstream bulk Lorentz factor. The ion 
spectrum (upper panel) is consistent with a thermal dis- 
tribution in the whole range of 70 we explore (from 70 = 3 
to 70 = 50), but for 70 < 5 (green for 70 = 5, yel- 
low for 70 = 3) the ion Maxwellian is much hotter than 
f° r 7o ^ 10 (black for 70 = 15, blue for 70 = 30, red 
for 70 = 50). Correspondingly, the electron distribution 
(lower panel) , which peaks at roughly the same energy as 
ions for 70 > 10, shifts to much lower energies for 70 < 5. 



Lyubarsky (20061 argues that the synchrotron maser 
instability should only operate for 70 > (m^/TOe) 1 / 3 . In 
this regime, the shock-compressed field grows on a scale 
smaller than the electron Larmor radius, and the non- 
adiabatic motion of incoming electrons into the shock will 
result in a ring-like distribution in phase space, which is 
synchrotron-maser unstable. It follows that a powerful 
electromagnetic precursor will only be generated in high- 
70 shocks. Our results show that, for rrii/m e ~ 16, the 
transition occurs at 70 ~ 10. For 70 > 10, when the 
precursor wave is strong, electrons reach energy equi par- 
tition with ions ahead of the shock. As discussed in §4.2| 
powerful electrostatic oscillations generated in the wake 
of the precursor dissipate their energy by heating a frac- 
tion of the incoming electrons, that eventually populate 
the high-energy comp one nt seen for 70 > 10 in the elec- 

Instead, for 70 < 10, wakcficld 



tron spectrum of Fig. 17 



oscillations are much weaker, so that the electron high- 
energy component tends to disappear (green for 70 = 5, 
yellow for 70 = 3). Also, since the precursor is less pow- 
erful in low-70 flows, ion-to-electron energy transfer in 
the upstream is suppressed, and the downstream electron 
spectrum shifts to lower energies. Yet, the downstream 
electron population still contains an appreciable fraction 
(~ 15%) of the initial ion energy, due to significant elec- 
tron energization by the cross-shock potential. 

A note of caution is required concerning the elec tron 
high-energy tail for 70 > 10. As shown in Fig. [XTJ this 

died 



component grows in time up to a maximum (reached be 
tween w p ji ~ 2000 and uj pi t ~ 3000), and then it recedes 
(although we cannot exclude a cyclical behavior). The 
growth time scales with the wavelength of wakefield oscil- 
lations, which is longer (in units of c/w p ;) for larger bulk 
Lorentz factors ( Lyubarsky 2006 ) . It follows that at fixed 



it = 2250, while 70 = 15 (black curve in Fig. [T7[ ) is al- 
ready at maximum, the electron high-energy population 
is still building up for 70 = 30 (blue curve) and 70 = 50 
(red curve) . At later times (not shown in Fig. |l7| ) , the 
electron spectrum for 70 = 30 and 70 = 50 will roughly 
resemble the case 70 = 15 presented here. 

In summary, our results for 70 = 15 are a good proxy 
for highly relativistic flows (70 > 10), but they do not 
apply to mildly relativistic shocks. 

6. SUMMARY AND DISCUSSION 

We have explored by means of 2.5D PIC simulations 
the internal structure and acceleration properties of rel- 
ativistic collisionless shocks propagating in a magnetized 
electron-ion plasma. This work complements the study 
of SS09, that analyzed relativistic shocks in magnetized 
electron-positron plasmas. To investigate how the shock 
properties depend on the conditions of the upstream flow, 
we vary the bulk Lorentz factor 70, the magnetization a 
and the angle 9 between the upstream magnetic field and 
the shock normal. Our results, confirmed by extensive 
convergence studies, can be summarized as follows: 

• For a ~ 0.1, the magnetic field is so strong that 
charged particles are constrained to move along the 
field lines, which are advected downstream from the 
shock. If the magnetic obliquity 9 < 6 CI n ~ 34° 
("subluminal" shocks), downstream particles can re- 
turn upstream by following the magnetic field. In this 
case, the electric current of "returning" ions can trig- 
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ger the so-called Bell's instability (Bell 20041, which 
generates powerful circularly-polarized Alfvenic-type 
waves in the upstream. The shock results from the 
nonlinear steepening of such waves, which happens 
quasi-periodically every few ion gyro-times (the so- 
called "shock reformation" process). Instead, for 9 > 
#crit ( "supcrluminal" shocks), no particle can propa- 
gate ahead of the shock along the field, and counter- 
streaming instabilities are suppressed. In this case, 
the shock is mediated by magnetic reflection of the in- 
coming flow of f the shock-compressed field (e.g., Alsop 



& Arons 1988J). In short, the dichotomy between sub 
luminal and superluminal configurations, presented by 
SS09 for electron-positron flows, still holds in electron- 
ion shocks. 



• If particle acceleration requires repeated bounces back 
and forth across the shock, as in the standard Fermi 
picture, only subluminal shocks should be efficient ac- 
celerators. Or, as measured in the upstream frame, the 
magnetic field needs to be within an angle ~ 6* cr it/7o — 
34° /70 from the shock normal. In subluminal shocks, 
particles are accelerated mostly via non-resonant inter- 
actions with Bell's waves (a Fermi-like diffusive pro- 
cess), although for 30° < 9 < 9 cl a acceleration by 
the upstream motional electric field may be important 
as well (Shock Drift/Surfing Acceleration; see SS09 for 
similar conclusions). In the downstream spectrum, the 
shock-accelerated particles populate a power-law non- 
thermal tail. For our representative case 9 = 15°, the 
ion tail contains ~ 5% of ions and ~ 30% of ion en- 
ergy, and its slope is —2.1 ± 0.1. The tail stretches 
linearly with time to higher Lorentz factors. Electrons 
enter the shock with a significant fraction (~ 15%) of 
the initial ion energy, yet their Larmor radius is still 
small compared to ions, and they are rapidly advected 
downstream with the magnetic field. It follows that 
the acceleration of electrons is less efficient (~ 2% by 
number and ~ 10% by energy, with a steep tail of slope 
—3.5 ± 0.1). In the upstream, the pressure of shock- 
accelerated ions (in fact, few electrons propagate back 
upstream) can significantly alter the structure of the 
shock, forming a smooth density "precursor" ahead of 
the actual "subshock" (this is the characteristic profile 
of "cosmic-ray modified" shocks). The overall width of 
the shock layer increases with time at ~ 0.1 c. 

• In sup erluminal shocks, the syn chrotron maser insta- 
bility (Hoshino & Arons 1991) generates a train of 
transverse electromagnetic "precusor" waves propa- 
gating into the upstream, that push on the incoming 
electrons, making them lag behind ions. This initiates 
longitudinal "wakefield" oscillations in the upstream, 
that transfer energy from ions to electrons ( |Lyubarsky| 
2006[ ). In high-obliquity shocks, electrons enter the 
shock with bulk energy comparable to ions. So, elec- 
trons acquire most of their energy in the upstream 
rather than at the shock, in contrast to what is as- 



sumed in mod els that only conside r the cross-shock 



potential (e.g., Gedalin et al. 2008). A few percent 



of the incoming electrons receive an extra amount of 
heating from the dissipation of wakefield oscillations. 
They eventually populate a prominent high-energy fea- 
ture in the downstream electron spectrum, which may 



resemble a (short) power-law tail. However, we remark 
that it r esults from hea ting, and not acceleration (as ar- 
gued by Hoshino|2~0 08 ) , and it does not extend in time 
to higher energies. Rather, as the incoming electrons 
get hotter, the gener ation of electromagnetic precursor 
waves is suppressed ( Hoshino fc Arons] [1991), and the 
electron high-energy component recedes to lower ener- 
gies. Hereafter, a cyclical behavior is likely, though 1111- 
demonstrated. Despite the significant wakefield heat- 
ing experienced in the upstream, which may facilitate 
injection into a diffusive acceleration process, no sig- 
nature of Fermi-accelerated electrons is seen in mag- 
netized superluminal shocks (see also SS09). 

• The properties of magnetized subluminal shocks do not 
significantly depend on the upstream bulk Lorentz fac- 
tor, in the relativistic regime 70 > 3 we explore. For 
superluminal angles, the sy nchrotron maser i nstability 
is suppressed for 70 < 10 (Lyubarsky 2006). So, for 
low-70 shocks the electromagnetic precursor is weaker, 
and transfer of energy from ions to electrons in the 
upstream is much less efficient. Yet, downstream elec- 
trons still account for as much as ~ 15% of the initial 
ion energy, due to significant electron energization by 
the cross-shock potential in the ion foot. 

• For lower magnetizations (a < 10 -3 ), the incoming 
plasma can filament, and the shock is mediated by 



the ion We ibel instability (Wcibcl 1959; |Medvedev fc 
LoebJ1999||Gruzinov & Waxman|l999|) , wRTcF ensures 



strong electron heating and efficie nt particle acceler- 
ation, for both ion s and electrons (Spitkovsky 2008a; 
Martins et al.|2009 1 . The nonthermal tail in the down- 
stream ion spectrum contains ~ 5% of ions and ~ 20% 
of ion energy, and its slope is —2.5 ±0.1. A similar 
tail is present in the electron spectrum, with compa- 
rable slope but a smaller particle (~ 2%) and energy 
(~ 10%) content. For magnetizations a < 10 -3 , the 
fields generated by the Weibel instability are stronger 
than the background field, and no difference persists 
between subluminal and superluminal configurations. 

• An important question is whether our 2.5D simulations 
can capture the relevant three-dimensional physics of 
electron-ion shocks. We performed a limited number 
of 3D experiments with relatively small computational 
grids (~ 3 ion skin depths along each transverse di- 
mension), obtaining essentially the same results as in 
our 2.5D simulations. Yet, large 3D runs are certainly 
desirable to confirm the overall picture presented here. 

Our findings may constrain the composition, Lorentz 
factor and magnetization of jets in Active Galactic Nu- 
clei (AGNs) and Gamma-Ray Bursts (GRBs)0 Syn- 
chrotron and inverse Compton emission from the core of 
blazar jets is usually attributed to high-energy electrons 
accelerated in mildly-relativistic internal shocks (70 ~ 2 
in the jet comoving frame). Such shocks should be quasi- 
perpendicular, since polarization measurements of radio 

19 Pulsar Winds are thought to be dominated by electron- 
positron pairs, and they will not be discussed here (we refer to 
SS09). If a significant ionic component is present, electron (and 
positron) acceleration may procee d via the Resonant Cy clot ron Ab- 
sorption mecha nism discussed bylHoshino et al. 1(19921 and Amatol 
|fc Arons|j2006) . 
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knots ( jGabuzda et al.|2004| |Pushkarev et al.||2005[ ) indi- 
cate that the magnetic field is mostly transverse to the 
jet axis. Electron acceleration in internal shocks is also 
invoked to explain the prompt emission of GRBs (e.g., 
Firan 12004] ) . If the jet flow carries a substantial toroidal 
held, internal shocks in GRBs will be magnetized and 
quasi-perpendicular. For both AGNs and GRBs, it is in- 
ferred that the emitting electrons should contain a sub- 
stantial fraction (e e ~ 10%) of the flow energy, and their 
spectrum is usually modeled as a power law extending 
over several decades in energy. 

We find that significant heating of electrons is a uni- 
versal property of relativistic magnetized shocks. The 
mechanism involved may vary depending on magneti- 
zation and magnetic obliquity, but post-shock electrons 
account for at least ~ 20% of the bulk pre-shock en- 
ergy, regardless of the upstream conditions. Although 
the agreement between this lower limit and the observa- 
tional requirement mentioned above is surely encourag- 
ing, we remark that most of our downstream electrons 
belong to the thermal component, whereas the emitting 
particles from AGNs and GRBs are usually thought to 
be nonthermal. Yet, our findings are important to con- 
strain the characteristic downstream electron energy to 
be > 0.2 jQiriiC 2 , with the following two implic ations 



First, this rules ou t models of GRB emission (e.g., |Bykov 
& Meszaros 1996[ ) which assume that only a small trac 



tion of electrons are accelerated up to high energies (and 
share the whole e e ), whereas the bulk of electrons remain 
cold (at ~ 7oTO e c 2 ). Second, the minimum Lorentz factor 
of a hypothetical power-law tail will be > 0.2'yomi/m e . 
For models of GRB and AGN emission that require a 
power-law distribution starting from smaller Lorentz fac- 
tors, this would suggest that electron-positron pairs may 
be a major component of the flow (see SS09). 

Regarding nonthermal acceleration of electrons, we 
find that quasi-perpendicular magnetized shocks are in- 
efficient particle accelerators, in contrast to what is re- 
quired in models of GRBs and AGN jets. However, 
the high-energy component that we find in the elec- 
tron spectrum of high-obliquity shocks, which actually 
results from heating (and not acceleration), may be 
(mis)interpreted as a short power-law tail. In fact, this 
component might reasonably explain the emission from 
AGN jets, which often requires an electron distribution 



extending o nly for one decade in energy ( Celotti fc Ghis-| 
ellini 2008). Otherwise, our findings imply that the jet 
pre-shock flow should be weakly magnetized (a < 10 -3 ), 
or maybe seeded with strong small-scale magnetic tur- 
bulence, that could facilitate the Fermi proce ss even in 
supcrluminal shocks ( Sironi fc Goodman|2007 ) . Alterna- 
tively, electron acceleration may occur not m shocks but 
in reconnection layers of a Poynting-dominated jet (e .g., 
Lyutikov fc Blandford||2003[ IGiannios fc Spruit]|2006[ ) . 

Finally, it internal shocks in AGlNs and GRB s are the 
sources of Hig h Energy Cosmic Ray ions (e.g., Sironi & 
Socrates 2010), then the pre-shock magnetic field needs 
to be nearly aligned with the shock normal, within an 
angle ~ # c rit/7o — 34° /7 (as measured in the upstream 
fluid frame) , for the shock to be "subluminal" . Alterna- 
tively, the acceleration of cosmic ray protons may occur 
at weakly magnetized shocks, like GRB external shocks 
or the termination shock ("hot spot") of AGN jets, or 
in the jet "sheath", where interaction of the jet with the 
surrounding medium may create subluminal configura- 
tions. In any case, we find that ion acceleration is always 
accompanied by substantial magnetic field amplification, 
especially ahead of the shock. In the presence of acceler- 
ated electrons, the upstream region may then contribute 
considerably to the observed emission. Synthetic spec- 
tra extracted from PIC simu lations of shocks (as done in 
unmagnetized pair shocks by Sironi fc Spitkovsky|2009b|) 
may possibly reveal characteristic signatures that could 
observationally distinguish such upstream emission from 
the downstream-based models that are usually employed. 
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APPENDIX 
A) DEPENDENCE ON THE MASS RATIO 

The reliability of our results may be questioned based on the artificially small mass ratio (mj/m e = 16) that we 
employ in most of our simulations. However, as we show here, for the parameters considered in the present work, a 
mass ratio of mj/m e = 16 is already large enough to separate the electron and ion scales, and capture the relevant 
physics of electron-ion shocks. 

Fig. [l8| compares the downstream particle spectra for three representative magnetic obliquities (9 = 15°, 30° and 
75°, witn 70 = 15 and a = 0.1) between two values of the mass ratio: mi/m e = 16 (dashed lines) and nii/m e = 100 
(solid lines). The spectra are computed at the same time (in units of w"; ) and at the same distance behind the shock 
(in units of c/u) p i). Very good agreement is obtained in all cases, for both ions and electrons. 

The downstream electron spectrum peaks at much higher energies than expected in the absence of any ion-to-electron 
energy transfer (compare with black vertical lines; dashed for mj/m e = 16, solid for mj/m e = 100), and the location of 
the peak is roughly independent of the mass ratio. It follows that the "effective" electron mass, including the average 
electron Lorentz factor, is a substantial fraction of Joirii, independent of nii/m e (for mj/m e > 16). The electron 
physics will then evolve on the ion time and length scales, which explains the similarity in Fig. H8] between electron 
spectra with different mj/m e , in terms of both the low-energy thermal part and the high-energy component. 
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Fig. 18. — Comparison of downstream particle spectra (upper panel for ions, lower panel for electrons) between different mass ratios: 
mi/m e = 16 (dashed lines) and m^/m,. = 100 (solid lines). The black vertical lines (dashed for m;/m e = 16, solid for m.j/m e = 100) 
correspond to ■yom e /mi, i.e., the location expected for the peak of the electron distribution in the absence of any ion-to-electron energy 
transfer. 



Various tests, for mass ratios up to uii/m e = 
7o), confirm the overall picture presented here. 



fOOO and for different upstream parameters (we changed both a and 



B) DEPENDENCE ON THE FIELD ORIENTATION WITH RESPECT TO THE SIMULATION PLANE 

We now assess the robustness of o ur results upon the orientation of the upstream magnetic field Bo with respect 
to the 2D simulation plane. In Fig. 19 (upper panel for ions, lower panel for electrons), we plot with red lines the 

ip = 0°), whereas black lines are for B lying in a plane 



In Fig. 

downstream spectrum if B is in the simulation plane (i.e., 
perpendicular to the simulation plane (i.e., tp = 90°). 

For low-obliquity shocks (e.g., 6 = 15° in panel (a)), our results are independent on the angle tp. However, for shocks 
in which the obliquity is close to the critical boundary 6 cl it — 34° between subluminal and superluminal configurations, 
we observe a substantial suppression of SDA for in-plane magnetic fields, both in electron-positron and electron-ion 
flows (see 8 = 30° in Fig. 21(b) of SS09, for pair shocks). As discussed in SS09, the case with tp = 90° is in better 
agreement with 3D simulations. For this reason, in the whole range of subluminal angles presented in the current 
work, we employ tp = 90°. 

For high-obliquity shocks (e.g., 6 — 75° in panel (b)), the agreement between in-plane and out-of-plane results is 
also remarkably good. The only appreciable difference is in the low-energy part of the ion spectrum (upper right 
panel). For tp = 90°, the plane perpendicular to the magnetic field (where downstream particle orbits would tend to 
lie) is almost degenerate with the simulation plane, resulting in a downstream ion distribution which is closer to a 2D 
(rather than 3D) Maxwellian. Instead, for tp = 0° the two planes are well separated, and ion isotropization in three 
dimensions is favored. For this reason, we preferentially employ tp = 0° for superluminal shocks. 

C) ION COUNTER-STREAMING INSTABILITIES 

In both unmagnetized and subluminal magnetized shocks, the interaction between the incoming flow and the shock- 
accelerated particles that propagate back upstream is a potential source of counter-streaming plasma instabilities. In 
this section, we determine how the nature of the instability depends on the upstream parameters (bulk Lorentz factor 
7o and magnetization a, choosing 9 = 0° as a representative subluminal obliquity). Primed quantities are measured 
in the upstream fluid frame, unprimed quantities in the downstream (simulation) frame. 

Let Ccr = n CR /rii be the ratio of cosmic-ray number density to number density of incoming ions. By "cosmic rays" , we 
mean the shock-accelerated ions moving ahead of the shock (we neglect the returning electrons, since their contribution 
is usually smaller). In the upstream frame, £ ( ' R = 7q Ccr; given that the population of returning ions is roughly isotropic 
in the simulation frame. Similarly, we define e CR = U CK / {^miriiC 2 ) as the ratio of cosmic ray energy density to kinetic 
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Fig. 19. — For fixed magnetic obliquity (9 = 15° in panel (a). 8 = 75° in panel (b)), downstream particle spectra at w p ;t = 2250 for 
different orientations of the upstream magnetic field Bo with respect to the simulation plane: magnetic field lying cither in the simulation 
plane (tp = 0°, in red) or in a plane perpendicular to the simulation plane (ip = 90°, in black). 



energy density of the injected ions, such that the average cosmic ray Lorentz fa ctor is 



= 7o(<WCch.)- In the 

upstream frame, e' CR — 7q£cr and 7^ R = e' CR /(' CR — Jo Jen- Informed by Figs. 14 and 16 we can assume to a first 
approximation that £ 0R and e CR (as measured in the downstream frame) do not significantly depend on 70 or a . 

The maximum growth rate of Bell's instability (wg dl ~ l/2(^ / CE uj p i) is achieved for a characteristic wavelength 
A' Bell ~ 47r(V(7 7 /CcR)c/wpi ( Reville et al. 2006 ). Here, cr' is the magnetization paiameter as defined in the upstream 



frame (cr' = 7q cr, for 9 = 0°). The corresponding Larmor frequency of background ions is oj 1 - — 70 w C i- 

Bell's instability governs the evolution of the shock if the following two conditions are simultaneously satisfied: (i) 
its characteristic wavelength Ag ell is smaller than r' L CR , the Larmor radius of cosmic ray ions; (ii) its growth rate is 
smaller than the Larmor frequency of background ions (wggn < w ci)i i- e -j the upstream plasma is magnetized. Since 
r L cr = 7cr/ c/uj p i, the first condition can be rewritten as 2 cr' < e' CR in the upstream fluid frame, which corresponds, 
in terms of quantities in the simulation frame, to 

2a<7o£ C R ■ (CI) 

When the above inequality is broken (either for large cr or small 70, at fixed e CR ), the physics of the instability 
changes. The dominant mode is now generated via gyro-frequency resonance with the cosmic rays (as opposed to 
Bell's instability, which produces non-resonant modes) . We observe this change of polarization in the dominant mode 
for magnetizations a > 0.3 (with fixed 70 = 15) and for bulk Lorentz factors 70 < 2 (with fixed a — 0.1). 

The condition (ii), that the upstream ions stay magnetized on the growth time of Bell's instability, can be rewritten 
as Ccr ^ m the upstream frame. In the simulation frame, 



7oCcr 



< 



For lower magnetizations, or higher bulk Lorentz factors, the upstream plasma can filament, 
instability dominates. For 70 = 15, this transition happens in our simulations around cr ~ 10~ 2 . 



(C2) 

and the ion Weibel 
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